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ABSTRACT 



An experimentcl investigation of the effect of variotion 
in flow angle on convection heat transfer and frictional 
resistance in a tube bank is undertaken for three sys tematicolly 
varied geometries. Data taken by previous investigators for 
a fourth geometry is incluaed in the analysis. 

Design-oriented curves and equations derived from analysis 
of the results are presented. The curves give the value of an 
inclination factor which may be used as an additional factor 
in the modified Grimison equation. 

Similar design oriented curves are presented which give an 
inclination factor which is in effect the ratio of friction 
factor at any angle to that for cross flow. Equations for 
the friction factor at ninety degrees are derived. 

Extensive use is made of electronic Digital computers for 
aata reduction ana analysis. 
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Program 

Listing 

A 

AK 

ANRE 

ANNU 

ARGFAC 

B 

BK 

BNRE 

BNNU 

C 



NOMENCLATURE AND DEFINITIONS 



Engineering Definition 

Symbol 

A The inclination factor equation inter- 

cept in the convection heat transfer 
regression analysis. 

K Flow coefficient for a sharp-edged thin 

plate orifice. 

Nj^^ Reynolds number chrough the thin plate 

orifice . 



Nj^jy Experimental Nusselt Number, Cross flow. 

F^ Arrangement factor in the moaified 

Grimison equation. 

B Coefficient of Reynolas Number in the 

convection heat transfer regression 
analysis . 

K Flo’w coefficient for sharp-edged tnin 

plate orifice calculated by data 
reduction program. 

N^g, Reynolas Number in the tube bank v;here 

transport data is a function of TF 

Ny Experimental Nusselt Number, any angle 

- Quotient of AK and BK taken so as to 

always re less than unity. 
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Program 

Listing 


Engineering 

Symbol 


Definition 


CNRE 


^RE 


Reynolds Number in the tube bank where 
transport data is a function of TFF, 


CNNU 


^NU 


Nusselt Number calculated by Colburn 
equation. 


CONST 


K 


The coefficient of Reynolds Number in 
the frictional resistance analysis. 


CPD 


c 

P 


Specific heat of dry air at constant 
pressure. 


CPH 


c 

p 


Specific heat of water at constant 
pressure. 


CPW 


c 

p 


Specific heat of wet air at constant 
pressure . 


DHF 


— 


Pressure drop across the thin plate 
orifice. 


DIAFAC 

DNRE 


iW° 


Two tenths power of the ratio of 
hyaraulic diameter to tube diameter. 

Floating point mode equivalent of NRE. 


DP 


— 


Pressure drop across the tube bank in 
inches of water. 


ENRE 


^RE 


Reynolds Numbers used in computation 
of friction factor from frictional 



resistance analysis 
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Program 

Listing 


Engineering 

Symbol 


Definition 


EXPO 


X 


Exponent Jfcf Reynolas number in the 
frictional resistance analysis 


FRICT 


f 


Friction factor computed from 
parameter developed by frictional 
resistance analysis. 


FRIFAC 


f 


Friction factor computed in data 
reduction program. 


- 




Convection heat transfer coefficient. 


IDO 


d 


Diameter of the thin plate orifice 


IRON, RNCD 


- 


Run code number for identification. 


Equation 


k 


Thermal conductivity of air 


KASES 




Maximum number of data sets in any 
run of) the data reduction or frictional 
resistance analysis programs, index 
of controlling DO Loop in each. 


LAF 


- 


Back pitch indicator. 


LANGLE 


- 


Angle indicator. 


LGF 




Arrangement indicator index of DO Loop 
representing the sequential number 
of a set of data being processed in 
data reduction program. 


M 


- 


Exponent of Reynolas Number in convection 
heat transfer analysis. 
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Progrcim 

Listing 

NRE 



NUMB 



Pi 



PB 

RHO 

SPHU 

TAVG 



TE 



TF 



Engineering Definition 

Symbol 

- Index of DO loop in frictional 
resistance analysis for computation of 
frictional resistance from parameters 
computed in the program. 

- Inaex analysis program of DO Loop 
in frictional resistance. Maximum 
number of data points in a aata set, 
in the frictional resistance analysis. 

Pj^ Pressure above atmospheric ahead of 

the thin plate orifice 

Pj^ Barometric pressure 

Air density 

X Specific humiaity 

Average or bulk temperature of air in 
the tube bonk. 

T2 Exhaust temperature; temperature of 

air after passing through the tube 
bank . 

Film temperature, mean between bulk 
temperature and tube surface 
temperature . 

Film temperature, for frictional 
resistance calculation, based on 
Tg - C (LMTD) 



TFF 
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Program 

Listing 

THUMB 

TT 

V 



WAO 



XCON 

Y 
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Engineering Definition 

Symbol 

- Floating point moue equivalent of 

NUMB 




log^K 

Y 



o 



Tube surface temperature 

Intermediate quantity calculated in 
the data reduction program. 

’.Veight rate of flow of air through 
the orifice. 

Natural logarithm of CONST. 

Net expansion factor, ratio of the 
flow coefficient of a gas to a 
liquid at the same Reynolds Number. 

Intermediate quantity calculated in 
the data reduction program. 

Dynamic viscosity of air. 

Angle of inclination of tuoe rank; 
cross flow is ninety degrees, 
parallel flow zero cegrees. 
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DEFINITIONS OF COMPUTER ORIENTED TERMINOLOGY 

Fixed-point variable - a variable written without a decimal 

point using the digits o, 1 9; 

may be positive or negative. 

Floating-point variable - a variable written with a decimal point 

using the digits o, 1....9, v;ith an 
optional proceeding plus or minus sign. 
The variable value may contain a 
multiplier v.hich is ten to some power, 
the multiplier being represented by 
the symbol E and the power of ten. 

DO loop - a method of performing the same 

calculation with different cata by 
repeatedly executing the statements 
following the DO statements the number 
of times inaicated by an inaex. 

BLOCK DIAGRAM SYMBOLS 
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Processing a group of program instruc- 
tions which perform a processing 



function 
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Input or output - any function of any 
input or output device. 



Decision or Logic - The Decision or logic 
function documents points in the program 
where a branch to alternate points is 
possible. 

Predefined process - a group of operations 
not detailed by the particular program. 



Terminal - the beginning end or a point 
of interruption of a program. 



An off-page connector 



Magnetic Tape 



Connector 
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RUN CODE NUMBERS 

A system of coding data points to reflect in a single 
multi-digit number, tube geometry, angle of inclination and 
sequential number of the test run was developea. The run code 
number is a four digit number; the first number indicating 
tube geometry, the second indicating tne angle of inclination, 
and the third and fourth numbers indicating the sequential 
number of the run. 



Tube geometry codes: 

Code Description 

1 Staggered tubes, ST = SL = 1.5 (Smith-Kiss and 

Bond -/v' allin data) 

2 In-line tubes ST = SL = 1.5 

3 Staggered tubes ST = 1.8 SL = 2.0 

4 In-line tubes ST = 1.5 D, SL = 2.0 



Angle codes: 



Code 

3 

4 
6 
7 
9 



Description 

30 degree angle of inclination 
45 degree angle of inclination 
60 degree angle of inclination 
75 degree angle of inclination 
90 aegree angle of inclination 



example: 1704 

tube geometry 1, ST = SL = 1.5 D^, staggered tubes, 
seventy five degree angle of inclination. 



fourth run 






4 
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INTRODUCTION 



While adequate empirically derived design information is 
available for heat transfer equipment where the tube banks are 
parallel to or normal to the flow of the fluid being heated or 
cooled, little information exists for the case where the tubes 
are inclined at some angle to the fluid flow. The case of 
tubes at an angle to fluid flow is not uncommon purticularly 
in marine boiler convection banks. It is the purpose of this 
thesis to investigate the effect of variation in flow angle on 
convection heat transfer and frictional resistance in a tube 
bank using a model scale apparatus, and to aevelop cans tc.nts 
v/hich will be useful to the designer of hi_at transfer equipment. 

The relationship originally developed by Grimiscn [1] * and 
later mouified by the transport data obtained by Justi and 
Leuder [2] v;ill be us d as the basic convection heat transfer 
relationship wh',^re the tube banks are inclined at ninety degrees 
tc thio tluid flow. A factor will be developed i.hich, when 
induced in the .moaificd version [3]of this c.,u.tici will give 
it the form: 






0,2 9/ P h f f 

' A- :r D ^t a 






1/3 



The itictor ie will oe obtained Irem aesJ. in curvcc, . nd 'wuen 

used will give a moaified value of Nusselt Number h.d, for 

t t 

k 

tubes inclined at an angle 0 to fluid flov;. 



so 

the 



Concurrently, with the investigation of tne effect of 
variation in flov/ angle on convection heat transfer tne effect on 
frictional resistance as expressed by the Fannihg type friction 
factor will be determined. 

•Numbers in brackets ur.signate reference^ lisLud ot the ena of text. 
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The availability of electronic computers makes possible 
more accurate reduction and analysis of experimental data 
than has been possible before. In addition, they permit a 
greater number of calculations to be carried out in a short 
time. Thus a systematic series of geometries will be tested 
by varying the longitudinal tube pitch and tube arrangement. 



Ill 
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THEORY 

Effect on Convection Heat Transfer 

The transfer of heat by a fluid passing through a tube 
bank is by the mechanism of forced convection in turdulent flow. 
The quantity of heat transferred through the stagnant fluid film 
or thermal boundary layer contiguous to the tube is a function 
of the following quantities: 

a) a characteristic length, taken in this development as 
the outside diameter of the tubes, D^, 

b) the coefficient of convection heat transfer, h^ , between 
the fluid and the tube surface. 

c) the thermal conductivity, k, of the fluid. 

d) the specific heat, c^, of the fluid at constant pressure. 

e) the dynamic viscosity of the fluid. 

f) the mass flow of the fluid, G, 

The convection heat transfer coefficient may be expressed 

i 

as a function of the other quantities: 

h = f (Dt, K, cp,yw , G) 

The choice of a characteristic length and of a temperature 
at which to evaluate those temperatures dependent physical 
properties is open to some discussion. When fluia flows across 
a single tube the outside diameter of the tube is the most 
logical choice, however in a tuoe bank the flow is effected by 
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the clearance between tubes and other definitions of character- 
istic length have been suggested. 

To preserve consistency with Piersonl" ^Gra...iioon^^ una Huge 
and because other definitions provide no proven advantage the 
outside diameter of the tube will be used. Similarly various 
definitions of the temperature at which to evaluate physical 
propecties have been suggested. 

That which is consistent with the work mentioned above 
and which is most logical for forced convection heat transfer 
is the temperature of the fluid film. As there is a variation 
in fluid film temperature through its thickness the film 
temperature is defined as the arithmetic mean between the bulk 
temperature and the tube surface temperature. Likewise the 
bulk temperature varies through the tube bank and is defined 
as the arithmetic mean of the inlet ana exhaust temperatures. 

The exact mechanism of forced convection heat transfer is 
not fully under sto'.d so that while analytical methods can 
sometimes be used, semi-empiricol or empirical relationships 
are more reliable. The oevelopment of empirical relationships 
is assisted by the derivation of dimensionless groups of the 
physical qualities discussed above. 

These groups may be derived by the methods of dimensional 
analysis. The applicable dimensions are: 

L for length 
T for time 
K for mass 
H for heat energy 
0 for temperature 



U > 
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As heat is equivalent to mecnanical v;crk it can be 

2 —2 

replaced by K force) x (resistance) or H • ML T . The 
dimensions of the remainder of the physical properties may be 
expressed as: 



D 



t 



h 



f 



L 



k 

c 

P 

/U 



MLT ^<f) 



-1 

-1 




G = 




-1 



If the Buckingham Pi Theorem is applied, the following matrix 
is developed: 



M 

L 

T 



0 

1 

0 

c 



1 

0 

-3 

-1 



K 

1 

1 

-3 

-1 



P 

0 

2 

-2 

-1 



/4 

1 

-1 

-1 

0 



1 

-2 

-1 

0 



From this matrix three dimensionless relations whicn nave 
proved useful in heat transfer work are aerived. These are: 



^f^t = L 

^ 1 

HT D 



= Nusselt number, 



which may pnysically be interpreted as the racio of the 
characteristic length to a stationary fluid layer conducting 
heat at the same rate under the same temperature aitference; 



I 

,1 

I 
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r^n » p — 1 

t = ML~ T~ L = Reynolas number, 

^ -1 -1 

ML T 

which may be interpreted as the ratio of inertia force to 
viscous force and: 



^ . ML~^T ^ = Prandtl Number 

which is the ratio of two constants of molecular transportation 
The relationship is usually written in terms of power functions 
as : 



hfD^ 



= A 




where the coefficient A and the exponents a and b are constants 



By correlation of empirical data where the fluid flow is 
normal to the tube bank the above equation has been shown to be 



hfDt 



= 0.292 





1/3 



where F^ and Fj^ are empirical correlations known as arrangement 
factor and depth factor respectively. The arrangement factor 
is a function of tube arrangei ient , staggered or in-line, 
longitudinal and transverse spacing and Reynolds Number. The 
depth factor is a function of the number of rows in the tube 
bank. The depth factor oetermined by Grimison becomes unity 
for ten or more rows. 
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tube 



Colburr!^^ has shown that for fluid flow parallel to the 
surface : 



K 



= 0.023 




1/3 



Where D„ is the hydraulic diameter. In oraer to be able to 
make a direct comparison with flow normal to the tubes the 
above expression may be written in terms of tube diameter. 



h^Dt 



0.023 






0.80 



D 

0 , 



0.20 



H 
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For cases where the fluid flow is neither normal nor 
parallel to the tube surface, but inclined at some angle theta 
to it, the rate of heat transfer as indicated by the value of 
Nusselt Number should lie somewhere between that for normal 
flow or cross flow, and parallel flow. It should be possible 
to express the effect of the inclined flow on forced convection 
heat transfer in one of two ways either as: 



1) a modification to the exponent of Reynolas Number, i.e.; 



hfDt 



=0.292 




0.60 f(©) 




where f 

2) 



(©) represents a function of the angle of inclination, 
as an additional multiplying factor: 



or 






hfDt 



0.292 F^Fj^ f (0) 




0.60 




1/3 
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An analytic preaiction of the effect of flow angle on 
convection heat transfer would be at best uncertain. As the 
earlier work with which this work is to be consistent is an 
empirical correlation, no theoretical prediction will be 
attempted. Instead, data derived from the test apparatus 
will be correlated witn the empirical equation to proauce a 
modifying factor which is a function of the angle of inclination. 



The frictional resistance of a fluid flowing across a 
tube bank is customarily treated by an equation f the general 
Fanning type wherein the pressure drop is treated as Joeing 
caused by flow through a series of orifices formed by successive 
major restrictions in the flow path. By methoas of dimensional 
analysis similar to those described above for convection heat 
transfer a dimensionless friction factor, f, may be aerived 
which is a function of Reynolos Number only: 



Where the coefficient C ana tne exponent X are constants. 
This equation may be modified to permit the use of convenient 
quantities without reducing its soundness: 



Effect on Frictional Resistance 



f = C 




f = 10.84 (10)^/^DP 
2 

NG 



) 

2 



I 
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where : 



10.84 = a constant which incluaes the effect of the 

acceleration of gravity in ft. per hr. per hr. 
and a head conversion factor from pounds per 
square foot to inches of water 



DP 



N 



= air density 

= pressure drop across the tube bank in inches of 
water 

= number of major identical restrictions. 

f in the case of a tube bunk N is equivalent to ; 



i the number of rows of tubes 



J 



- mass fiow 



The question of what temperature should be used to aetermine 
the properties of the fluia again arises. The temperature has 
a significant effect on gas density and viscosity which directly 
effect the friction factor and tne corresponding Reynolas Number. 
Utilization of the mean bulk temperature in the bank has been 
shown to be unsatisfactory by previous investigators. Empirical 
aevelopments led to the establishment of the convention that 
the physical properties be evaluated at a temperature equal 
to the tube surface temperature minus eight tenths the log mean 
temperature difference for staggered arrangements or equal to 
tube surface tem.perature minus nine tenths the log mean 
temperature difference for in-line arrangements. The minus 
sign applies where heat is transferred from h- ated tubes to 



a cooler fluid 



■ 
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As is true with the convection heat transfer tne 
prediction of variation of friction factor with flow angle 
does not lend itself to analytic prediction. Accoraingly no 
attempt at analytic prediction will be made. 

The choice of a purely empirical approach while not as 
intellectually stimulating as a theoretical development has 
the distinct advantage of producing results which are 
reproduceable and can be used with some confidence in the design 
of heat transfer equipment. 
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EXPERIMENTAL TECHNIQUE 

The technique followed in the investigation is that of 
forcing air through a bank of heated tubes which are arranged 
in various geometries and inclined at various angles to the 
air flow in order to measure changes in parameters indicative 
of the quantities of heat transferred ana the frictional 
resistance in the bank. Data is obtained for three geometries, 
all with a pitch normal to the air flow of 1.50 tube aiameters. 
Pitches in the direction of flow are 1.50 and 2.00 tube diameters 
and the tubes are arranged in staggered and in-line configura- 
tions. For each geometry, measurements are taken with the 
tube bank at angles of 90, 60, 45 and 30 degrees to the air 
.flow. The angle is varied by inclining the tube bank in a 
test section of constant cross section and length. 

The tubes are internally heated by steam, the inlet pressure 
being controlled so that the tube surface temperature is very 
near 212 degrees Farenheit. Ambient air is heated by being 
forced through the tube bank by a steam turbine driven blower. 
Variations in the mass flow of air are obtained by varying the 
speed of the supply blower. The mass flow of air is measured 
by means of a square edged orifice and appropriate manometers. 
Temperature sensors in the air duct upstream and downstream 
from the tube bank are used to aetermine inlet and outlet 
temperatures. Tube surface temperatures are determined by 
sensors located on the tube surfaces tnroughout the bank, 
generally on the upstream sector of the tube. Tube surface 







i 
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4 
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temperatures used in the computations are then the arithmetic 
mean of the measured tube temperatures. The principle parameter 
for determining the effect of frictional resistance is the 
pressure drop across the tube bank measured by manometers. 
Transport data is evaluated at appropriate temperatures with a 
humidity correction obtained by standard psychometric techniques 
applied to the specific heat. 

In addition to the data obtained from the three geometries 
tested, the data determined by previous investigators has also 
been included. All this data has been placed through the 
computation process and analyzed for correlation of both the 
previous work and the extended investigation. 



26 



PROCEDURE FOR Rfc:DUCING LXPERIHLKTAL DATA 
The purpose of the reduction of experimental data is the 
calculation of Nusselt Numbers, Reynolds Numbers and friction 
factor. This process may be divided into four parts: 

1) calculation of weight rate of flow through the apparatus 

2) calculation of Nusselt Number and the corresponding 
Reynolds Number 

3) calculation of friction factor and the corresponding 
Reynolds Number 

4) calculation of Nusselt Numbers by the modified Crimison 
equation ana the Colburn equation. 

This section develops in detail the calculation procedure, 

1) Weight rate of flow: 

The weight rate of flow calculation is perhaps the most 
complex of those performed. The equation is: 



WAO = weight rate of flow 
AK = flow coefficient 
Y = expansion factor 
g = gravitational constant 
P = air densit^^ 

( P 1 -P 2 ) = pressure differential across the orifice 



The method used to obtain weight rate of i low from orifice 
dimensions and pressure differential is a standard engineering 
technique and h^s been developed in detail by both Bond ano 



V.AO = A^ (AK) 




where: 
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By combining constants and applying proper unit conversion 
factors the equation reduces to: 

WAO = o. 02416 (IDO)^ (AK) Y / DHF 



where : 

IDO = orifice meter diameter in inches 

DHF = pressure drop across the orifice meter in inches of oil 



The expansion factor Y is a function of diameter ratio, 
ratio of differential pressure to absolute static pressure and 
the ratio of specific heats. As shown in Appenaix D the factor 
Y cdn be calculated for each data point by applying an equation 
of the form 

Y = 1 - (DHF) 

p^ - 13.596PB 

The flow coefficient, AK , is a function of orifice diameter, 
uiameter ratio, and pipe Reynolds Number. The aerivation of 
the expressions for flow coefficient in Appendix D gives 
equations of the form: 

AK = C 2 + 

✓ ANRE 

However Reynolus Numoer in the pipe is given by the 
expression : 

ANRE = y (IDO ) (AK) C /dHF 
where : 



ANRE = Pipe Reynolds Number 
V = Expansion Factor 



IDO 



Orifice Diameter 
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yU = Dyneimic Viscosity for air, the expression for which 
was developed in Appendix D and is * 

^ = (.16 TI + 112.6) xlO~^ 

AK = flow coefficient 
C = a constant 

DHF = differential pressure across the thin plate orifice 
meter. 



It can be seen that flow coefficient and Reynolds Number 
are mutually dependent and an iterative process must be used 
to determine satisfactory values for these two quantities. 
This can be done by first calculating intermediate quantities 
composed of those factors other than AK which appear in the 
Reynolds Number expression specifically: 



V 



2.239 

0.16TI + 112. ( 



X 10^ 



and : 

y » 1 (dhf) 

p^ + 13.596PB 

A trial value of AK depending upon orifice diameter may be 
assumed and a corresponding value of Reynolas Number may be 
calculated : 

ANRE = V X Y X AK 

This value of Reynolds Number may be used to calculate a 
second value of flow coefficient which in turn may be used to 
compute a new Reynolds Number. This process can continue 
until the comparison of the new value of AK with tne previously 
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calculi. ted value shows that the increase in accuracy doi_s not 
warrant repetition. This final value of flow coefficient 
can then be used to colculaue weight rate of flow of air 
through the apparatus. 

mo = (86.978) (IDO)^ (AK) (X) 



When the weight rate of flow has been determined certain 
other values may be determined. These include: the average 

temperature in the tube bank in aegrees Rankine: 

TAVG . TI. i _TE ^ 

and the specific heat of wet air by the method aetailed in 
Appendix D, the log mean temperature difference: 



XLMTD = ~ TI 




ana film temperature: 



TF 




TT 



2.0 

When these quantities have been computed everything 
necessary for the second major calculation step, the computation 
of Nusselt Number has been completea. The computation of 
Nusselt Number may then be undertaken. 



2) Nusselt Number: 

Nusselt Number may be defined as 

Q 



Nusselt Number 



S (XLMTD) 



Dt 

K 
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where : 

Q = total heat transferred 
Dt = tube diameter 
S = total heating surface 
XLi-.TD = log meun temperature difference 
k = thermal conductivity of air 
The relation for Nusselt Number may then be expressed as: 

ANNU = NAO (CPf) (TE - TI ) 

LTT(XLMTD) (0.000024TF + 0.01344) 

where : 

^ = vvAO (CPv.) (TE-TI) 
k= 0.000024TF + 0.01344 
Dt/S = effective length L x TT 

by utilizing the interior dimension of the box, minimum possiole 

tube length, and combining constants the equation reduces to: 

ANNU = ..AO (CPW) (TE - TI) 

XLMTD (0.01158TF + 6.4829). 

For flow in which the tubes are inclined the surface area is 
increased. The surface area at any angle of inclination is 
equal to the surface area at ninety degrees divided oy the 
sine of the angle of inclination, thus for angles other than 
ninety degrees the expression for Nusselt Number becomes: 

BNNU = WAO (CPW) (TE - TI ) SINE 0 



XLMTD (0.01158TF + 6.4829) 



ANNU (SINE ») 
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This is the proper value of Nusselt Number for inclined flow 
and is the one result of the second part of the calculation. 

It is now necessary to calculate the Reynolds Number through 
the test section, Reynolds Number may be expressed as: 

BNRE = WAO (Dt) , 
where : 

A = free flow area through the test section, a projected 

area which remains constant though the angle is varied 

— 7 

= dynamic viscosity of air = (,16TF x 112,6) x 10 

The equation reduces to: 

BNRE = WAO X 10^ 

(0,046256TF + 32,552) , 

With these two results, Nusselt Number and Reynolds Number the 
second part of the computation is complete, 

3) Friction factor 

For the frictional resistance calculation some quantities 
which have not been previously calculated are required, 

2 new definition of film temperature is necessary. This 
definition is a function of tube temperature and log mean 
temperature difference and may be expressed in the form: 

TFF = TT - Cg (XLMTD) 

where is either 0,800 or 0,900 depending upon whether the 
tube geometry is staggered or in-line, A Reynolds Number for 




\ 
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use with the friction factor may be calculated using the 

N 

value of film temperature just computed: 

CNRE = WAO X 10^ 

(0.046256 TFF + 32.552), 

The friction factor is expressed as: 

FRIFAC = 10.84 (RHO) DP 10 ^ 

N^ VJAO^ 

The density” of air is calculated using the equation: 

p 

— — = RT in the form 
-^RHO 

RHO = + 13.596PB) 5.1815 / 

53.349 TFF = 459.69 

When RHO has been obtained the friction factor may be calculated 

FRIFAC = 1.084 (RHO) (DP) 10 ^ 

WAO^ 

0.1673^ 

4) Nusselt Number by the Grimison Equation 

When the friction factor has been obtained, the 
computation continues to calculate the Nusselt Number by the 
modified Grimison equation. 

The Grimison equation is expressed as: 

TNNU = 0,292 (ARGFAC) F^ (BNRE) ^ CPW/Y j^^^ 

As the test section is ten rows deep the depth factor, 

F^ , is unity and can be eliminated from the expression. The 
quantities Reynolds Number, BNRE, and specific heat, CPW are 
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available from previous work, and analytic expressions in 
terms of temperature exist for dynamic viscosity and thermal 
conductivity. It has been shown in Appendix D that the 
arrangement factor may be expressed as a function of Reynolds 
Number for specific geometries. With proper suosti tutions the 
expression for Nusselt Number becomes: 



TNNU 



0.292 ARGFAC (BNRE) 



(0.16TF +112.6) (10 



.C00024TF + 0.01344 



.j'hen this value of Nusselt Number has been determinec, a 
factor FTHLTA is defined by forming the quotient of the 
experimental Nusselt Numoer BNNU, and the Grimison Nusselt 
Number , TNNU : 

FTHLTA = BNNU 
TNNu. 



When FTHLTA has been computed, the computation of Nusselt 
Number by the Colburn expression is undertaken. The similarity 
in the Colburn expression and that for tne Grimison equation 
leads to simplifying trie calculation by calculating the Colburn 
Nusselt Number using the previously culculated Grimison 
Nusselt Number: 

|/3 

where DIAFaC is a conversion factor tne tune aiarneLer to 
hydraulic diameter. 



CNNU = 



,02 3 (BNRE)^*“^(DIAFAC) I CP 



0.60 



TNNU 



292 (BNRE) 



( aRGFAC ) 



K 



73 










I 



I 

I 

i l' 

I l' 

< 




I 



I 



t 



34 



Dividing CNNU by TNNU ana solving for CNNU: 

CNNU = , 0 7877 (Br,RE) ^ ( DIAfAC ) TNNU , 

ARGr AC 

This form is readily digestible by the computer, /.'hen ChNU 
has been computed the calculation of one dota point is 
complete , 

'.Vhen the calculations have been completed for all data points 
the processing of experimental data is complete and analysis 
can commence. 
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ANALYSIS OF CONVECTION H£AT TRANSFEH RESULTS 

Figures 1 thru 4 are composite plots showing the 
experimental variation of Nusselt Number with Reynolas Number 
for each of the four geometries which have been tested. The 
variation in slope of the plots as the angle of inclination 
changes indicates that the inclination factor is not a 
simple multiplying constant but is an exponential function of 
Reynolds Number. 

Fg is defined as the ratio of the experimentally determined 
Nusselt Number to the Theoretical Nusselt Number, calculated 
by the modified Grimison equation, at a particular value of 
Reynolds Number. Thus for every experimental point a value 
of Fq as Reynolas Number may be plotted. 

Rather than plot these values by hand, the data processing 
facilities at David Taylor Model Basin were utilized to enhance 
the speed and accuracy of the process. An existing program 
was modified and used to provide the needed output. This 
program is described in David Taylor Model Basin Report No. 

2037 fyj and in Appendix G. 

The program takes as input the values F^ and the corres- 
ponding Reynolds Number plus identifying codes which serve to 
individualize each of the points. The output is the equation 
of three lines of regression which are produced by means of 
the method of least squares. In aadition, there is a graphic 
output which plots the three lines of regression which are 
determined. Samples of these output forms are included in 
Appendix G. A complete summary of the regression analysis 
coefficients is shown on page 42. 
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COMPOSITE PLOT 
V/^Ri£~^lON OP NUSSELT number 
WITH REYNOLDS NUMBER 
s^=-s,^i.5Dj sro. 






Ujr 



REYNOLDS NUMBER X 10 



F! 6 . / 
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COMPOSITE PLOT 
VARIATION OF NUSSELT NUMBER 
WITH REYNOLDS NUMBER 

S^= S = l.5Dr IN-LINE 




FI 0. 2 
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COMPOSITE PLOTI 
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Sl= 2 . 00 t 
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Of the three regression equations for each geometry, one 
represents consideration of all the points, one is determined 
by all the points within two standard errors, ana the third 
represents the line determined by points within one standard 
error, where a standard error for a set of N points is given 
by the relation: "■ ^ 



The choice of which equation to utilize was made through 
inspection of the SC4020 plots of the regression equations. 

In almost all cases, the values of at low Reynolds Numbers 
were outside the two standard error criteria, and, since the 
first three points are among those used to make a first deter- 
mination of exponent, these points influence the phase of the 
regression line unduly. Further, the data for values of Reynolds 
Number less than 2000 is suspected to be influenced by laminar 
flow. Thus only that data which fits within the two standard 
error criteria was allowed to determine the regression equation 
and the second set of coefficients were used for further analysis. 

For design work, the equations as shown in Table I may be 
used if the designer is working at a specific angle of inclination 
and geometry. Entering with Reynolds Number as an argument, 
a value of F^ will result which can be used as an additional 
factor in the Modified Grimison equation. 







where Y avg. is 
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For greater convenience, two aesign plots have been 
preparea. Figures 9 ana 10 in tne cnclusion section, give 
the factor F^ as it varies witn angle at tnrec constant values 
of Reynolds Number. There are 12 curves, three for each of 
the four tube arrangement geometries analyzed. 

Entering with the angle of inclination ana the geometry, 
the designer may read a value of F^ for each of the two values 
of Reynolds Number which bracket the argument he requires. 

By linear interpol tion, F^ can be particularized to the 
desired Reynolds Number. 
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TABLE I 



Summary of Inclination Factor Regression Coefficients 
For Points Aithin Two Standard Errors 

Fg = A * B 







A 


B 


M 


S^=1.5D^ 


45 


0.5087726 


1.056159(10)"^ 


0.6510998 


Sj^=1.5D^ 


60 


0.6856272 


4.110530(10)"'^ 


0.9889 


STG 


75 


-0.6487573 


0.6318375 


0.1069 




90 


0.8435844 


1.62334(10)"'^ 


0.7754198 


S^=1.5D^ 


30 


0.6208538 


2.209353(10)"^ 


0.9889 


S^=1.5D^ 


45 


0.7091857 


3.031199(10)"^ 


0.9889 


IN-LINE 


60 


0.7860375 


4.310665(10)"^ 


0.9889 




90 


1.155682 


-4.277594(10)"^ 


1.00999 


S^=1.5D^ 


30 


7.290824(10) 


4.853093(10)"^ 


0.2980699 


S^=2.0Dt 


45 


0.7260181 


3.172470(10)"^ 


0.9889 


STG 


60 


0.7262158 


1.74172 3(10)"'^ 


0.8012598 




90 


0.7653439 


7.028655(10)"^ 


0.4712899 


S^=1.5D^ 


30 


0.4970187 


3.129761(10)"^ 


0.9889 


Sj^=2. OD^ 


45 


0.6867613 


2.621991(10)"^ 


0.9889 


IN-LINE 


60 


0.6011853 


3.337236(10)"^ 


0.9889 




90 


0.8751028 


2.385977(10)"^ 


0.9889 
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ANALYSIS OF FRICTIONAL RESISTANCE 

Figures 51 through 66, included in Appendix H, are plots 
of friction factor against Reynolds Number for indiviaual 
geometries and inclinations. Several points pertinent to 
subsequent analysis are observable from these curves. The 
non-linear behavior of the full logarithmic plot at low Reynolds 
Numbers is apparent in these plots. This behavior is a 
characteristic of small tube apparatuses that is consistent with 
the observations of previous investigators and is probably 
caused by laminar flow in the tube bank. In order to keep the 
analysis as general as possible data points with Reynolds 
Numbers less than 2500 have been omitted. Additionally it can 
be seen from Figure 55 that the data for the thirty degree 
angle, = 1.5D^, In-Line, is completely inconsistent with 

all the remaining data and is therefore omitted from the 
analysis. The inconsistenty is undoubtedly due to faulty 
instrumentation or operator error during the collection of 
this data. 

The friction factor can be expressed analytically as a 
function of Reynolds Number by fitting the data to an equation 
of the form: 

FRIFAC = CONST (CNRE)^^^*^ 
using the method of least squares. 

A computer program, described in Appenaix H, has been 
utilized to compute these parameters which describe the friction 
factor as a function of Reynolas Number for a particular geometry 
and angle of inclination. The equations which include these 
parameters are summarized in Table II. 
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TABLE II 

FYPO 

f = CONST (CNRE) 





G 


CONST 


S,j,-S^-1.5D^ 


90 


1.931187 


STAGGERED 


75 


0.445318 




60 


0.045465 




45 


0.65543 


S^=Sj^=1.5D^ 


90 


0.261858 


In-LINE 


60 


0.037337 




45 


0.100758 


S^=1.5D^ 


90 


0.431623 


Sj^=2.0D^ 


60 


0.187797 


STAGGERED 


45 


0.070171 




30 


0.040784 


S^=1.5D^ 


90 


0.297399 


Sj^=2 . OD^ 


60 


0.074656 


IN-LINE 


45 


0.017001 




30 


0.033258 



EXPO 

-0.457575 

-0.045756 

-0.196590 

-0.043053 

-0.081586 

-0.119514 

-0.024007 

-0.086865 

-0.028274 

-0.089713 

-0.040784 

-0.055798 

-0.05331 

0.210122 



0.060968 
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Figures 5 thru 8 are composite plots showing the variation 
in friction factor with Reynolds number at different flow 
angles for particular geometries. The plots are based on the 
equations listed above. From these plots it may be seen that 
the set for on 3^ = 5^^ = 1.5D^, STAGGERED, (Fig. 5) is not 
consistent with the remaining data. The data for this geometry 
was taken by previous investigators v;ho experienced some 
difficulty in correlating it. This previous data will be 
omitted from further analysis. 

From the equations representing the variation of friction 
factor with Reynolas number, it can be seen that there is no 
apparent pattern in the variation of either the coefficient or 
exponent with angle of inclination for particular geometries. 

Friction factors for Reynolds Numbers of three thousand 
through eignt thcusana, in increments of one thousand, were 
computed as part of the analysis. Program Values of Reynolds 
Number = 3000, 5000 and 8000 are tabulated below to indicate 
the variation with flow angle for particular geometries. 

REYNOLDS NUMBER = 3000 

Friction Factor 

GEOMETRY 30 45 60 75 90 

St = Sl = 1.5D^, STG - 0.0464 0.2194 0.a087 0.1140 
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REYNOLDS NUMBER = 3000 continued 







Friction Factor 




GEOMETRY 


30 


45 


60 


75 


90 


S,p=S^=2 • OD^ , J.L. 


- 


0.0831 


0.0972 


- 


0.1363 


Sy=1.5D^, S^=2.0D^,STG 


0.0648 


0.1312 


0.1498 


- 


0.2153 


S^=1.5D^, Sj^=2.0D^, I.L. 


0.0542 


0.0914 


0.1144 


- 


0.1902 


REYNOLDS NUMBER = 5000 








S^=S^=1.5D^,STG 


- 


0.0454 


0.2426 


O.3016 


0.1095 


St=Sl=i»5DtjI.l. 


- 


0.0821 


0.1033 


- 


0.1307 


S^=1.5Dy, S^=2.0D^ STG 


0.0668 


0.1258 


0.1476 


- 


0.2060 


S^=1.5D^, S^=2.0D^I.L. 


0.0599 


0.1018 


0.1176 


- 


0.1850 


REYNOLDS NUMBER = 8000 








S^=S^=1.5D^ STG 


- 


0.0445 


0.2661 


0.2952 


0.1054 


^T=Sl=1*5D^ I.L. 


- 


0.0812 


0.1093 


- 


0.1258 


S^=1.5D^, S^=2.0D^ STG 


0.0686 


0.1206 


0.1457 


- 


0.1977 


S^=1.5D^, S^=2.0Dy I.L. 


0.0575 


0.1124 


0.1206 


- 


0.1801 


From these values, 


a relationship 


can be 


developed wni 



expresses friction factor at any angle as a function of tne 

friction factor at ninety degrees by simply taking ratios of 

the values and plotting them: i.e. 

R© = fQ 
f9© 

where f© = the friction factor at any angle, 

f9© = the friction factor at ninety degrees as expressed by 
f, = CONST (CNRE)^^^° CONST and EXPO being 
particularized to the geometry. 
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The values so derived are tabulated below: 



= 3000 
RE 

GEOMETRY 




F© 

30 


45 


60 


90 


Sy=S^= 1.5D^ I.L. 




- 


0.6097 


0.7131 


1.000 


5^=1. 5D^ S^=2,0Dip 


STG 


0.3010 


0.6093 


0.6958 


1.000 


Sy=1.5D^ S^=2.0D^ 


I.L. 


0.3500 


0.4805 


0.6015 


1.000 


Nre = 5000 


1.5D^ I.L. 




- 


0.6282 


0.7904 


1.000 


S^=1.5D^ S^=2.0D^ 


STG 


0.3243 


0.6107 


0.7165 


1.000 


S^=1.5D^ S^=2.0D^ 


I.L. 


0.3388 


0.5502 


0.6356 


1.000 


"re = ®°oo 


S^=S^= 1.5Dy I.L. 




- 


0.6455 


0.8688 


1.000 


S^=1.5D^ S^=2.0D^ 


STG 


0.3470 


0.6100 


0.7370 


1.000 


S^=1.5D^ = 2.0D 


■£ I.L. 


0.3182 


0.6241 


0.6696 


1.000 



Values of the ratio F0 are plotted and are incluaed in the 
Conclusion section as Figures 11 through 13. By entering the 
applicable curve with values of flow angle the ratio of 
friction factor at the angle € to that at ninety degrees may 
be obtained. For values of Reynolds Numbers other than those 
tabulated, linear interpolation may be used. 
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CONCLUSION 

Convection Heat Transfer 

Figures 9 and 10 represent the plotted results of the 
equations developed by the regression analysis of individual 
inclination factors. These plots show the variation of 
inclination factor with angle for in-line and staggered 
geometries. Six curves are plotted in each case representing 
constant values for Reynolas Number of 2000, 5000 and 8000 and 
particular values of transverse and longituainal tube pitch. 

The curves can be used directly as a design aid by applying 
the value of inclination factor obtained as a multiplier to 
the modified Grimison equation. Values of the inclination 
factor for Reynolds Numbers other than those plotted may be 
obtained by linear interpolation. The modified Grimison equation 
now becomes: 

It should be noted that the variation of inclination factor 
with angle of inclination is by no means linear so that there 
are regions for both staggered and in-line geometries in which 
relatively large changes in angle of inclination cause small 
changes in the rate of heat transfer. The reverse situation 
also can be obtained. The choice of which region is optimum 
must be left to the user who is familiar with the requirements 
and limitations of his particular job. 
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FRICTIONAL RhSISTANCE 

Figures 11 through 13 represent the plotted results 
of the analysis of frictional resistance for each of the 
three geometries for which reliable data was obtained. Each 
figure contains three curves, one each for values of Reynolds 
Numbers of 3000, 5000 and 8000, Additionally on each figure 
there appears an equation in the form: 




This is the equation for the friction factor in each case. 

Two uses of the curve are suggested. The first is to 
calculate an absolute value of friction factor. This method 
should be used with caution for while the method of obtaining, 
reducing and analyzing frictional resistance data is believed 
to be more accurate that that used by previous investigators 
it does not agree closely with their work. The second use 
is to estimate qualitatively the effect of varying the flow 
angle on the frictional resistance in a particular piece of 
heat transfer equipment. 
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RECOMMENDATIONS FOR FURTHER INVESTIGATIONS 

The apparatus, test procedure and method of reducing 
and analyzing data in this thesis represent a suitable method 
of empirically determining the effect of flow angle on 
convection heat transfer and frictional resistance in a tube 
bank. The results presented are limited in that they include 
only four possible tube geometries. 

It is recommended that the series be extended to induce 
all values of transverse and longitudinal tube pitches likely 
to be found in modern heat exchangers. Such work would involve 
manufacture of new tube banks and sealing strips. It is 
furth> r recommended that data oe taken at an angle of seventy- 
five degrees both for other geometric. s which may be tested and 
for those for which this angle has not been previously used. 

The additional point wcula be of great assistance in firmly 
locating the upper end of the inclination factor curves. 

If a blower of larger capacity should become available 
the test data for all geometries should be extended to higher 
Reynolds numbers. Prospective investigators are advisea that 
the limit of the currently installed pressure instrumentation 
has been reached by the present air supply system. They 
are also aavised that the manufacturing processes are long 
and tedious and that tr.e condition of the laboratory steam 
plant should be carefully investigated before unoertaking to 
use it. 

The reduction ana analysis of test data could be facilitated 
by combining the tnree c>,,mputer programs used in this work into 
one program to reduce and analyze dt^ta in one computation. 
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APPENDIX A 
Apparatus 

At the onset of tne project, it was felt that the apparatus 
as designed by Bond ana Wallin and later used by Kiss and Smith 
woula not be adequate for the task without considerable 
modification. The apparatus as constructed ana availadle haa 
been usea for one basic geometry by Kiss ana Smith and had 
an inherent disadvantage in that it was difficult to make any 
alteration in tube geometry. 

The design of a new test section was therefore undertaken. 
The principles which guided the design are: 

1, The section must be of the same cross sectional 
dimensions as the original apparatus as verified by 
Bond and VJallin. 

2, The section must have the capability of rapid ture 
geometry changes, 

3, All manufacturing operations necessary for completion 
of the design must be feasible with equipment avail- 
able in the Webb shop facility, 

4, The design must utilize the same tube banks that were 
prepared for the previous investigators. 

A design which satisfied the four stated requirements was 
developed. The configuration is generally the same as that 
used in the earlier experime^nts. The steam supply and drain 
line were used without moaif ication, as was the air supply 
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system including the orifice meter. The entire test section 
was rebuilt in accordance with the plan shown as Fig. 15 in 
this appendix. 

The test section consists of two rectangular wooden boxes, 
closed at one end, which serve as the entrance and exhaust 
portions of the test unit. Connecting these two end pieces to 
form a continuous square duct, is a center section of aluminum 
alloy. This center section is designed with stationary side 
plates and movable top and bottom plates. These features 
enables the tube bank tc be adjusted to any angle from vertical 
through 60 degrees from the vertical. This was the practical 
extent of adjustment in that the tube banks were not long 
enough between headers to permit another range. This is, 
however, considered to be a sufficiently wide range, as it 
permits angles of flow incidence from 30 degrees to 90 degrees, 
at any practical increment. This design parameter sized the 
length of the aluminum center section at 30 inches. 

The end sections are made up of 3/4 inch exterior grade 
marine plywood. The joints are made up with glue and joined 
with #8 X 1% " wood screws. The sides of the section are 
routed 3/16" by 3" to receive the sides of the aluminum center 
section such that the center plate will fit flush. The 

wooden sides are also routed at the top and bottom edge, 3/16" 

wide by I 4 " deep to receive the edges of the sliding top plates. 

The top and bottom of the vyooden intake ana exhaust sections 
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are routed to a depth of 3/16 inch to effect a flush fit for 
the sliding plates of the aluminum center section* They are 
also routed across the width and to a length which will 
accommodate the sliding plates through the range of inclination 
of the tube bank. This range is shown in Figure 17, the 

o 

drawing of the section. The intake section as a 6 inch 90 

✓ 

elbow secured to the end for delivery of air through a 6" 
diameter hole. The exhaust section has a similar stove pipe 
elbow made up to the side of the section. The heated air 
exhausts through a 6 inch diameter hole in the box then through 
the stove pipe exhaust line to the atmosphere. 

Both intake and exhaust sections were sanded smooth and 
finished with two coats of shellac prior to assembly. After 
assembly, all edges were caulked with a silicone rubber sealant 
and the stove pipe connections were made up with asbestos cloth 
gasketing using Permatex compound for a leakproof joint. 

The center portion of the test section is fabricated from 
2024-T3 aluminum alloy. This alloy was chosen for its free 
machining characteristics, since several machining operations 
were required in manufacturing the numerous sealing strips 
between the tube banks. It was later learned that the temper 
of this plate was not correct in that when the plates were cut 
to size and the edges and sliding grooves milled, the plates 
warped out of plane due to locked in rolling stresses. After 
several unsuccessful attempts to straighten the plates by 
different methods, the decision was maae to anneal the plates 
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the plates and thus stress r*. lieve them. Several industrial 
contacts were made, ana through the courtesy and kindness 
of the Sylvania-Corni ng plant in Hicksville, N.Y., the plates 
were annealed. Proviaed thus with plane surfaces, construction 
of the center section may proceed. One further note could be 
added for the benefit of those who may be curious, a more 
proper alloy choice for this application is 2024— TS3 this 
temper oeing stress relieved by stretching. The siae plates 
were made up with all edges milled and fastened to the wooden 
end pieces by flat heat wood screws, courtersunk to effect 
a plane surface. The top movable plates were maue up having 
all edges milled and two 3/16" wiae by 1/16 inch deep groove 
milled in the sides to effect a seal with the siae plates and 
to provide a channel for guidance of the moving plates. Each 
top and bottom plate (four total) has a projection along each 
edge v;hich slides in a channel. 

The channel was maae from an aluminum storm window shape 
and proved to have insufficient strength to hole the weight of 
the tube bank and secure the plates in position as well. To 
correct this, heavy steel yokes were made from angle and 
threaded to close the section ana secure the sealing strips in 
place. All but one of the channels were removed, the one 
remaining only to serve as a positioning device. Manufacture 
of the sealing strips represented the largest expenditure of 
time for any one operation. For each of the 12 combinations 
of angle and geometry, 22 sealing strips were required, 11 
each for top and bottom. Each sealing strip required edge 
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milling, groove milling, and 22 holes arilled or milled to 
close tolerances for completion. In all, 5680 holes were 
machined, and a total of 264 strips were preparea. The sealing 
strips proved very satisfactory in maintaining the tube geometry 
while providing excellent sealing characteristics. The only 
modification required was to provide a small hole in the end 
of each strip to permit the separate pieces to be wired to 
hole them in place while changing tube geometry. An isometric 
sketch showing the arrangement of the center test section 
and the sealing strips is provided as Figure 15. Drawings 
showing the sealing strip as a part of the test section and in 
detail appear as Figures 17 and 18, A photograph of the 
seeling strips in detail is available as Figure 22. 

The section is supported from a laboratory table by two 
wooden frames and a framework of steel angle v;hich was modified 
only to facilitate the changing of the geometry. This 
modification consisted of welding an angle on each corner of the 
metal structure to provide support for the tube bank headers 
when changing the sealing strips. 

The tube bank is thut used by previous investigators. 

It is made up of 10 rows of <5 inch O.D. copper tubes, 22 tubes 
per row, with each row of tubes being fitted with an inlet 
and exhaust header of Ih inch copper pipe, A detail drawing 
of an individual tube bank is shown in Figure 16. The inlet 
headers are connected to tl"ie steam supply moisture trap by 
3/8 inch copper tubing. Many of these lengths of tubing had 






{ 



I 




I 



68 



had to be replaced to permit the freedom of movement required 
in the tune bank. The exhaust headers were connected in the 
same manner to a drain header of inch copper pipe. An 
additional problem was that these lines kinked when moved aoout 
due to strain hardening and several were replaced for this 
reason. For some geometry changes, the exnaust ana supply 
lines were completely disconnected to avoid this problem. 

Two small holes are provided in one side of the center 
test section for thermistor leads. Four holes are provided 
for pressure taps in the test section. Two are in the inlet 
section, positioned in the middle of the vertical height and 
4 inches from the center section. Two taps are provided in 
the exhaust section on a similar vertical position, 4 inches 
from the test section. Thermistor probes are provided in the 
6 inch inlet piping just before the entrance elbow and in the 
6 inch outlet section to avoid radiation effects from the 
tube bank. 

Previous investigators had used various methods of sealing 
the banks; References 4 and 8* The method used in this series 
of tests, though long and painstaking in preparation, proved 
worth the effort while testing. There was no need to use any 
sealing compound once the tubes were in place and the top 
and bottom moving plates maoe tight. The tolerance to which 
the sealing strip holes were prepared was 0.008 inch and there 
was little or no leakage. Further, a geometry change could 
be made in about tv;o hours on the average and there was no glue, 
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sealer, etc. to clean from the tube to cause further delay. 
In this regard, the extra effort was rewaraed and this 
arrangement is recommended for later investigators. 

Figure 14 shows an overall isometric view of the steam 
plant and test apparatus. 
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DETAIL 'a' 

TUBE PENETRATION - SEALING STRIPS 
TUBES SHOWN AT 90“ AND 30°, Sl=2.0D. 



TOP PLATE 




DETAIL 



'B' - JOINTURE OF PLATES 
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Instrumentation 

Temperature 

Instrumentation for ottaining temperature readings was by 
thermistors installed in air intake and exhaust lines ana on 
tubes to sample temperatures throughout the tune bank. 

Thermistors are solid state electronic semi-conductors 
which exhibit a large negative temperature coefficient of 
resistance. Their small physical size and the large negative 
temperature coefficient of resistance make thermistors 
particularly well suited for application in small scale exper- 
imental models as units can be located within a very limited 
space and great accuracy can be obtained with a relatively 
simple electrical circuit. 

The thermistors used were Fenwal GB35J1 bead type thermistors. 
These units are 0.043 inches in diameter and have a nominal 

o 

resistance of 5K ohms + 20% at 25 C, The thermistors were 
calibrated by the use of a constant temperature bath. Boiling 
water and boiling benzene were used to obtain calioration 
beyond the ranged tiie bath. A seunple calibration curve is 
included as Figure 23 at tne end of this section. Thermistors 
used for inlet and outlet temperatures were fitted into 
probes made from brass tuning and sealed with silicone rubber 
sealant. The probes were located in the center of' the inlet 
and exhaust ducts at a point around a ninety degree bend from 
the tube bank so that they would be screened from direct 
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radidtion from the tube bank. Some difficulty was experienced 
in affixing the thermistor units tightly to the tute bank so 
that a true tube temperature would be inaicoted. The problem 
was finally overcome by taping the thermistors to the tune 
with plastic electricians tape. Prior to placing the 
thermistor on the tube, a small groove was ground on the tube 
surface to receive the bead. The tape and pigtail leaas were 
then covered with a high densa.ty epoxy. This effectively 
held the thermistor in place even during the curing period. 
Thermistor surface sensor assemblies consisting of thermistor 
beaus bonded to teflon or silicone impregnated fiberglass tape 
are commercially availaole. Tne cost of these sensor assemblies 
is approximately eight tim s that of tne thermistor units 
so that it was not considered economically feasible to use them. 
For similar reasons it was decided to undertake calibration 
of the units rather than to purchase unius which track a 
standard curve within narrow limits. 

The thermistor leads were connected to a selector switch 
which switched individual units into a Wheatstone Bridge circuit. 
The resistance necessary to balance the bridge was measured 
and this value was used to enter the calibration curves to 
determine temperatures in degrees Farenheit. A diagrammatic 
sketch of thf: temperature measuring circuit appears as 
Figure 24. 

A high failure rate was experienced with the thermistor 
units. Tne highest percentage of failures were due to 
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accidental severing of the small diameter wire pigtails at the 
point where the wire entered the glass beaa. These failures 
were caused simply by the rough service to which the units 
were subjected, particularly those installed on the tube 
surface. In addition, several thermistors were found to be 
erratic in behavior, culminating in a shiort circuit within 
the bead after a short time. No explanation for this performance 
has been found. Faulty units were replaced with spares until 
the supply was exhausted. At this time, it beccime necessary 
to replace the exhaust temperature thermistor with a mercury- 
in-glass thermometer. 

As an aside, it should be added that the thermistors were 
not given fair trial in this application. Tne necessity to 
move the tube bank about roughly placed undue strain on these 
units. Based on this experience, haray commercially prepared 
mounts such as those mentioned above are recommended. 
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Pressure 

Instrumentation for pressure was by manometers containing 
oil or water as appropriate. 

The manometer used to determine the pressure on tne 
upstream side of the orifice was an open water manometer. The 
barometer reading was converted from inches of mercury to 
inches of water and added to the manometer reading to obtain 
an absolute pressure upstream of the orifice. 

An oil filled U-tube difierential manometer was used to 
indicate the pressure drop across the orifice plate. A water 
filled U-tube differential manometer and a slant tube manometer 
were used to measure the pressure drop across the tube bank. 

A Bourdon tube pressure gauge and a bellows type compound 
gauge were used to measure steam inlet ana drain pressures 
respectively. 

Flow 

Instrumentation for weight rate of flow was by thin plate 
orifice meter with pressure tops located at one pipe diameter 
preceeding and one half pipe diameter following the inlet face 
of the orifice plates. To obtain/a large spectrum of Reynolds 
Numbers, three separate orifices of 2 inch, 3 inch, and 4 
inch diameters were used. The weight rate of flow was obtained 
from the differential head across the orifice. It was found 
that it was possible to cover the entire range of weight rates 
of flow which could be obtained with the available steam 
turbo-blower with the 3 inch orifice. 
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APPENDIX B 

Test Equipment Operating Proceoure 

Operation of the test equipment can be divided into two 
categories: 1) Lining up, ligliting off and operation of the 

steam plant and auxiliary equipment and 2) Operation of the 
test apparatus itself. 

The following procedure is followed in putting the steam 
plant into operation: 

1) The cooling tower is checked ^nd filled as necessary 
to ensure a proper quantity of cooling water. 

2) Cooling tower fans and the circulating pump are 
started . 

3) The condensate pump is started. 

4) Fires are lighted in the boiler following the procedure 
given in the boiler operation and maintenance manual. 

5) When the steam pressure has reached about 100 psi trie 
air ejectors are lined up and sufficient nozzles cut 
in to maintain about ten inches of vacuum in the 
condenser. The most satisfactory operation was obtained 
when both second stage nozzles and no first stage 
nozzles were used, 

6) The turbo blower is lighted off using a throttle pressure 
of about 5 psi. The blower lubricating oil tank should 
be checked to ensure that an aaequate supply of 
lubricating oil is flowing. 
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7) The steam separator is drained ana the tube bank 

steam supply and drain valves are openeu. A pressure 
of ten psi at the steam separator was found to be 
very satisfactory. 

When the above operations have been completed ana the 
conditions in the apparatus as reflected by inlet, exhaust, 
and tube temperatures have stdiLized the procedure enters its 
second phase, the collection of data. 

Test procedure consists primarily of varying the quantity 
of amounts of air supplied to the apparatus and recording the 
data provided for by the installed instrumentation, as aescrioed 
in Appendix A, 

The weight rate of flow of air to the test section is varied 
by throttling the blower supply steam with fine control being 
accomplished by manipuliition of the dischcirge bleed valve. 
Satisfactory operation is obtained by starting at low flow rotes 
and increasing gradually to liigh rate:^. v'-ith tnis procedure, 
the apparatus stabilizes quickly so that the interval oetween 
runs is very short. It should be noted that because the test 
section and the turbo blower are supplied by steam from tne 
same line fluctuations in boiler pressure caused by variations 
in blower demand will cause fluctuations in the pressure to the 
tube banks. As the tube temperatures ana exhaust temperatures 
are sensitive to the supply pressure variutions should be 
minimized by ensuring the tube bank supply is the same before 



each run 
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The pressure drop across the orifice is indicative of 
flow rate, accordingly successive runs may be maue by adding 
small increments of pressure drop by throttling the blower. 

During testing, it is advantageous to make a running 

plot of the quantity TT-TI \J DUF ver sus \J DHF on full 

TT-TE ! 

logarithmic paper. TT-TI \J DHF may be considered 

TT-TE 

directly proportional to Nusselt Number whiley' DHF may be 
considered directly proportional to Reynolds Number. A linear 
plot on full logarithmic paper then serves as a check on the 
consistency of the data. Care should be taken not to take 
such a plot as indicative ol the validity or accuracy of the 
data as the proportionality is not exact. A sample plot appears 
as Figure 25 at the end of this Appendix. 

Orifice plates of three different diameters are availoble 
for use with the apparatus. Utilization of the three inch 
orifice permits data to be taken over the entire range of 
flow obtainable with ti.e turbo blower. The other two orifice 
sizes were not used. 

If an immediate change of geometry is anticipatea upon 
completion of uata taking, allov;ing the blower to run after 
the steam supply to the tube bank has been secured cools the 
metal portions of the apparatus quickly and facilit(utes handling. 

Securing procedure is essentially trie reverse of the 
lighting off procedure. Upon securing in cold weather the 
cooling tower should be drained to prevent freezing of the pipes. 
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RUNS 2601-2626 




iTdhf 



FIG 25 
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APPENDIX C 

Development of Equations Used in Qfomputations 

I 

Computerization of the calculation procedure required the 
development of equations for a number of quantities which are 
usually tabulated or taken from curves in engineering 
applications. Development of equations which are applicable 
in the range of temperatures ana pressures which exist in the 
apparatus reduce the complexity of the program and permit a 
great saving in core storage within the machine. The quantities 
for which equations are developed herein are: 

1) Thermal conductivity of air 

2) Dynamic viscosity of air 

3) Specific heat of wet air 

4) Flow coefficients for orifice meters 

5) Expansion factors for orifice meters 

6) Grimison Arrangement factors 

7) Nusselt Number by Colburn equation 

The symbols used in this appendix are those developed for 
use in the computer program. 



1) The thermal conductivity of air varies in an essentially 
linear manner in the range of temperatures which occur in the 
apparatus. An equation may be written in the form: 

^ ^ X T = C 

Vv'here 

k' = thermal conductivity of air“^in Btu 

hr-ft °F 

i\lc = change in thermal cc nauctivity per degree Farenheit, 

,6t 
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T = Temperature in degrees Farenheit 
C = a constant 

Values of thermal conductivity used to develop the equation were 
taken from Table 2 of Gas Tables by Keenan ana Keys. [^9] Tne 
equation for thermal conductivity is: 
k .000024 (T-40) + .0144 
or 

k = 0.000024T + 0.01344 

This equation gives values whicn correspond exactly with the 
tabulated values betv^een 40 and 240 degrees Farenheit. 



2) Dynamic viscosity, like thermal ccnauctivi ty , vurics 
linearly with temperature in the range of application and can be 
expressed in the same form: 



= aa 

AT 

Vv'here : 



X T + C 



- dynamic viscosity in lbs 

sec-f t 



At 



= change in dynamic viscosity per degree Farenneit 



T = Temperature in degrees Farenheit 
C = constant 

Values used in the development of the equation were taken from 
Table 2 of Gas Tabl s by Keenan and Keyes. The equation for 
dynamic v;i ^cosity is: 



- j^» 



16 (T-40) + llij 
o>‘ 



.16T + 112.. 



] 



X 10 



X 10 

-7 



-7 



■ii 



m 
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3) The development of equations for specific heat of 

wet air at constant pressure requires the selection and adoption 
of applicable equations for the specific heats of air ana water 
and combining them with proper weigntlng. 

The equations selected were taken from Table I of 
Thermodynamics by Faires (l(^. 

For dry air: 

CPD = 0.219 + 0.342T x lO""^ - .293T x lO"® 

and for water: 

CPH = 1.102 - 33. IT'"^ + 416 

T 

where: T = average temperature in the tube bank in degrees 

Rankine tne specific heat of wet air is then determined by: 

CPW = SPHU (CPH) + 7000 CPD 
SPHU + 7000 

where: CPlv = specific heat of wet air 
CPD = specific heat of ary air 
CPH = specific heat of water 
SPHU = specific humiaity in groins per pound 

4) Flow Coefficients 

For square edged thin plate orifice flow meters with taps 
one pipe diameter upstream and one half pipe diameter doivnstream 
the following equation for flow coefficient obtains: 

AK = Ko + b 
where : 

AK = flow coefficient corresponding to any particular 
set of the parameters: pip'^ diameter (D), ratio of orifice 

diameter to pipe diameter 0,, and pipe Reynolds number (ANRl) 









I 
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Ko = Flow coefficient when pipe Reynolds number becomes 
infinite 

b = an empirical multiplier 
= 1000 
\/ ANRE 

Ko and b may be expressed in the follov;ing forms: 

Ko = (0.6014 - 0.01352 + (0.3760 + 0.07257 

0.00025 

+ .0025D + + 1.5^^^ 

b = (0.0002 + 0.0011 ) + (0.0038 + 0.0004 ) X 

D D 

1^^ + (16.5 + 5D^^^ )J 

For situations such as exist in the test apparatus where 
there are a limited number of possible values of pipe diameter 
and diameter ratios specific values of Ko and b may be 
calculated and applied to the basic equation. It should be 

1 s 

noted that although jQ represents a number less than unity 
raised to a relatively large power it cannot be ignored if 
acceptable accuracy is to be obtained. 

The values for hCo and b for the various orifice plates 
available are: 

Two inch orifice: 

Ko = 0.59781 
b = .00080 
Three inch orifice: 
ko = 0.6180 
b = 0.00133 
Four inch orifice: 

Ko = 0.6735 
b = 0.00229 
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The expressions for flow coefficients for each orifice then become 

Two inch orifice: 

AK = 0.59781 + 0.80 

ANRE 

Three inch orifice: 

AK = 0.6180 4 1.33 ^ 

,ANRE 

Four inch orifice: 

/ 

AK = 0.6735 + 2.29 
,'anre 

5) Expansion Factor 

The equation for the expansion factor through a square 

edged thin plate orifice meter wita taps one diameter upstream 

and one half pipe diameter downstream may be expressed as: 

V = 1 - (0.41 4 0.35/ ^ 

K 

where: /j! = ratio of orifice diameter to pipe diameter 

X = ratio of differential pressure to absolute static 
inlet pressure 

K = ratio of specific heats. 

■A 

If a constant value of K = 1.4 is assumed ana values of/J 
are calculated for the two, three and four inch orifice meters 
the equation reduces to: 



For 


2" orifice: 


Y 


= 1 


- .295X 




For 


3" orifice: 


Y 


= 1 


- .307X 




For 


4” orifice: 


Y 


= 1 


- .339X 




For 


utilization 


in 


the 


program the factor 


X must be 


expressed 


, in the same 


units 


as the input data. 


The numerator 



X is equal to the differential pressure across the orifice 

✓ 
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meter (DHF) which is in inches of oil. The denominator is 
equal to the sum of the barometric pressure (PB) expressed 
in inches of mercury and the pressure above barometric ahead 
of the orifice meter expressed in inches of water. All 

units will be converted to inches of water so the expression 
for X becomes: 

X = 0.827DHF 

+ 13.596PB 

where: DHF, P^ , and PB are input quantities as defined in the 
preceeding paragraph. 

0.827 = the specific gravity of the oil used in the manometer 
13.596 = the conversion factor from inches of mercury 
to inches of water. 

Inserting the expression for X and simplifying the expressions 
for expansion factor as used in the calculation are obtained: 

For 2” orifice: 

Y = 1 - 0.24464 DHF 

Pj_ + 13.596PB 

For 3" orifice: 

Y = 1 - 0.25457 DHF 

+ 13.596PB 

For 4" orifice: 

Y = 1 - . 28131 DHF 

+ 13.596PB 

6) Grimison Arrangement Factors 

The arrangement factors used in the Grimison Equation 
are a function of longitudinal ano transverse pitches, Reynolds 
Number and the arrangement (staggered or in-line) of the tube 
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bank, and are not amenable to a simple analytic expression. 
However as only one transverse pitch, 1.5 tube diameters, and 
a limited number of combinations of back pitch and arrangement 
were tested it was found feasible to plot arrangement factor 
as a function of Reynolas Number only for each geometry. The 
curves so plotted were parabolic in form. An equation of the 
form (BNRE)^ + C 2 (BNRE) + = ARGFAC could be vjritten for 

each of them by utilizing the LaGrange Method. The equations 
so developed are of the form: 



ARGFAC = (BNRE)^ + (BNRE) + 
where ; 



in 

GEOMETRY x (10)^^ x (10) 











IN-LINE 


.463 


-.2963 


1.063 


ST = S, = 1.5 D. 

J.J t 








STAGGERED 


-1.29 


.3462 


.94059 


St °t’ 


1.5 






IN-LINE 


5.7407 


-1.774 


1.10518 


= 2.0 D^, S^ = 


1.5 D^ 






STAGGERED 


4.074 


-1.574 


1.1698 



Nusselt Number by the Colburn Equation 



In order to establish a lov;er limit wnich would be comparaole 
to parallel flow the Colburn equation v;as used to calculate 
Nusselt Number. The Colburn equation is: 

K J 



K 



.023 



GD 



h 



A 



.8 



where : 



GD^ = Reynolds Number based on hydraulic diameter. 
7u~ 

For the Colburn Nusselt Number to be oirectly comparable with 
these developed experimentally or by the Grimison equation it 



I 
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must be based on tube diameter rather than hyaraulic diameter. 
To accomplish this both sides of the equation are raultij^lied by 
the ratio of tuoe ui. -meter to hydraulic diameter. 




in the nomenclature used 



CNIv’U = 0.023 (BNRE)^*®*^/ ^ 



D 



~-p¥ 



1/3 



H 



The hyaraulic diumeter is dofi.ned as four times tne free flow 
area diviaed by the wetted perimeter. For a s tagger eo tube 
arrangement ; 



Dh = 4 






TT4) 



n D, 



.J 












Dh = 4 



^T^L ^t 
T7 



Ol 



\ 



\ -i-- 



For the in-line arrangement: 



Dr = 4 



1 



L 



Vl - ^ 

TT D, 



D 



T 



J 



= 4 -(TO 



T L 



TTD, 






' 7' 






0 - 
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Since only one value of and two values of are considered, 
the hydraulic diameter can be numerically evaluated for these 
values. For = 1,5D^, = 1,5D^ 



= 4 (l.SDjj,)^ " °t^ 



DtTT 



D„ = 9- TT 



H 



D, 



TT 



Similarly for St = 1.5D^, D1 = 2,0D^ 



Dh 



The factor then can be evaluated for a particular 

K. 

geometry. This factor has been given the name DIAFAC and is 

included in the Colburn equation to give it the form: 



CNNU = .023 (BNRE)°*®° (DIAFAC) ^ CPW Aj 



a/3 
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APPENDIX D 

Experimental Data Reduction 
Program Description 

The computer program descriped herein is a FORTRAN IV 
program designed for use with the IBM 7090 data processing system. 
The program calculates Nusselt Number, Reynolas Number and 
frictional resistance coefficient fer varying flow angles from 
experimentally derived data in the Haeberle Laboratory, ..ebb 
Institute. The constants ana to some degree the method used are 
pc-.r ticularized to the apparatus used for data collection. 

The first step in the program is to calculate the sines of 
the angles of Inclination at which tests were conducted and to 
assign to them floating point variable names in order that they 
may be repeateoly called without recalculation of the sine 
function. 

Control of the program is achieved througn a DO loop 
utilizing the variable M which represents the sequential number 
of a set of aata being processed. It is required that the 
maximum value of M called KASES be supt lied as input each time 
the program is used. The value of M's incremented by one each 
time a set of data is processed. The DO statement refers the 
machine to a READ statement causing a new set of aata to be 
read in and processed. When the maximum value of M (M= KASES) 
is reached the program is completed. 

Required input to the program is: 

1) Run code number, a four digit fixed point varianle, 

( IRCN) 



I 
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2) Pressure arop across the thin plate orifice meter 

in inches of oil a three or four aigit floating point 
variable, (DHF). 

3) Pressure above atmospheric ahead of the orifice meter 
in inches of water, a three or four digit floating 
point variable, (PI). 

4) Barometric pressure in inches of mercury, a four aigit 
floating point variable, (PB), 

5) Pressure drop across the tube bank in inches of water, 
a three or four digit floating point variable, (DP). 

6) Temperature of the inlet air in degrees Farenheit, a 
three digit floating point variable, ( TI ) . 

7) Temperature of the exhaust air in aegrees Farenheit, 
a four digit floating point variable, (T£). 

8) Temperature of the tube surfaces in degrees Farenheit, 
a four digit floating point variable, (TT). 

9) The diameter of the orifice plate in inches, a single 
digit fixed point variable, (IDO). 

10) The specific humidity of the inlet air in grains per 
pound, a three digit floating point variable, (SPHU). 

11) A geometry factor indicating whether the tube geometry 
is staggered or in-line, a single digit fixed point 
variable, (LGF). 

12) An angle factor inaicating the angle of inclination of 
the tube bank, a single digit fixed point variable, 
(LANGLE) . 

13) Back pitch factor, indicating whether the back pitch 

is 1 . 5D or 2.0D, single digit fixed point variable (LAF) 
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The program first calculai.es the intermediate quantity 
V for use in the pipe Reynolus Number calculation. 

The program then queries the value of orifice diameter to 
determine the path to be taken in calculating the Reynolds 
Number thru the orifice. The constants chosen for calculating 
the expansion factor and flow coefficients are functions of the 
diameter of the orifice meter. The procedure followed is 
the same for all orifice diameters. A value of flow coefficient 
(AK) is assumed, and a value of V is calculated. 

The three quantities which have been calculated are 
multiplied together to determine an initial estimate of Reynolds 
Number . 

When an estimate of Reynolds Number has been made a better 
estimate of flow coefficient is obtained. 

The estimate of flow coefficient is then compared to the 
initial value. If the two values vary by one half of one per 
cent or less the calculated value of BK is taken as the Value 
of the flow coefficient and the value of Reynolds Number is 
accepted and the calculation continues. If the difference 
between the two values of flow coefficients is greater than one 
half of one per cent the calculated value of flow coefficient 
is used for an initial value for recalculation of Reynolds 
Number and flow coefficient. This process continues until the 
one half of one per cent criterion is reached. 

When satisfactory values of Reynolds Number and flow 
coefficient have been established the program proceeds by 
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squaring the orifice diameter and converting it to a floating 
point variable for use in the weight rate of flov; calculation, 
Weignt rate of flow is then calculated. 

The average temperature in the tube bank is calculated 
and converted to degrees Rankine. With this average temperature 
the specific heats of dry air and water are calculated and 
combined to determine the specific heat of the air in the 
apparatus , 

Log mean temperature difference and film temperature are 
then calculated. 

All quantities necessary for the calculation of Nusselt 
Number are nl!>w available and line quantity is calculated. 

The angle factor (LAf.’GLE) is used to select the proper value 
of the sine of the angle of inclination previously calculated. 

This sine function is then multiplied by the calculated Nusselt 
Number to determine the Nusselt Number for the angled configuration 

The Reynolus Number in the tube bank is calculated. 

The Program then queries the geometry factor to determine 
if the geometry is staggered or in-line and branches accordingly 
to calculate film temperature for the friction factor 
computation. 

Reynolds Number for use in the friction factor calculation 
is calculated. 

The density cf the air in the apparatus and the friction 



factor are calculated. 
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After the friction factor has been computed the program 
again queries the geometry factor and then queries tne back 
pitch factor to determine the geometry so that tne proper 
arrangement factor and diameter factor for use in the Grimison 
and Colburn equations can be computed. When the arrangement 
factor hcs been selectea anu comL uted it is employed in the 
Grimison equation to calculate a value of Nusselt Number. The 
value of Nusselt Number calculated by the Grimison equation is 
diviaed into that obtained by use of the experimental uata. 
value of Nusselt Number is calculated using the Colburn equation. 

All computations have been completed so that an output 
tape may be written. Output consists of: 

1) Run Code Numoer as reaa in (IRCN) 

2) Nusselt Number determined from experimental aata (BNNU) 

3) Nusselt Number determined from the Grimison equation (T NU) 

4) Nusselt Number determined from the Colburn equation (CNNU) 

5) The ratio between the two Nusselt Numbers (FTHETA) 

6) Reynolds Number based on the film temperature used in 
the Nusselt Number calculation (BNRE) 

7) Friction factor (FRIFAC) 

8) Reynolds Number using the film temperature used in tne 
friction factor calculation. 

Upon completion of writing the output data another set oi 
data is read in anu the program repeated. With tne 7090 system 
calculation of each set of data .equires about one seventh of 



one second 
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A block diagram of the program is presentea herewith. The 
block diagram is a graphic representation of the data processing 
system procedures used in the solution of the problem. The 
diagram is included to assist the reader with a visual record 
of operauional and logic steps in the program and the relaiionship 
each step in the program has to the solution. 

The symbols used in the aiagram are the standard IBM symbols 
representing machine functions. A Key to these symbols is 
included in the Nomenclature section. 

A program listing also follows the block diagram. 

The computer facilities used for these computations were 
those at the David Taylor Model Basin. Data was prepared using 
equipment available at Webb Institute. 
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IF ( N T ) I . * I F ,1 / 

V; K I T r . ( b » r ’ ) 

. 4 

si: AC ( 5 , ) I RCN , Dl- F • P i , , DM did d T d 3U , »PHU, LGF , L AnGL'C ,LAF 

V- ( s - ^ 3 s/ ( ( d 6 d I ) d 1 2 . 6 ) ) d I . 7 ) 

IF ( IGG-- ) 3d. 7 
A K = . F ' 

Y- ( I . ( . 2 44G4 *• CHF ) / ( P 1 r 1 3 . . ) ) < :>yj\n ( UHF ) 

AnFu.^ AK < V *FY 

v'- . 5'-^7r: I f .r ' /SORT ( A\ME ) 

Ot T 
AK= . F 2 

Y = ( 1 . - ( . 26467 K FHF ) / ( P I F 1 3 . H ’h ) ) < SOM T ( OMr ) 

Af\R> t=AK > V-5^ Y 

6X= . C I . ; F 1 . 3 3/ wPMT ( Ai-iFD ) 

CC F = ' 

AK-~ .68 

Y=.( 1 . . 28 1 31 « : HF )/ ( P 1 f 1 ) ) d .P I (OH' ) 

ANR ^ A’*' K V » Y 

CK= . 6 7 - f 2 . -/SGR r ( ANP‘ ) 

IF( AK-Erv )P, 1 d i 
C=PK/ AK 
CO fP P- 
C:= AH 

IF( .9PS-C) 1 d d3 
AK=, F 

IF < 1 I^O- 3 ) i 4 d r * P 

AK= .< 

tV d 1 f3C * 2 

SAld r . S I ^ AV « Y 

F AV G^- ( T I > TE ) / 2 • F 4 6 ^ - 6 9 

COC- .-liF( .?4^>dAVC/l. 4)~( . .'OS* ] A VC,4r _/l . M 

CPh - 1 . 1 ( 3 3. I /Sold { I A V ^ ) ) f ( ^ I . /TAV J ) 

CP>.-- ( ( S\ HUFCF H) 4(7.. . FCPO ) ) / ( "PHU d . . ) 

XLi'^ri;=( { -ri )/AL;d( ( rr-i i ) ) 

r F = ( ( T E 1 r i ) / .r . : F I M / : . V 

ANNld ( AU ^CPV: * ( fh - T I ) ) / ( XL, y 1 0 ^ ( • 1 d8 TP' F6 . 4 ^2 8 ) ) 

IF ( LANGur --4 ) 24 . 8 , 2< 

BNNU-= 5 IN3 



(- 'VC. r 



>^F-N 



SOb:'fCF STAT‘D vr. -4 r 



IFF ( S) 



3 / 



CL > r J 

25 f NNC; = Ar4.\iU * S 1 i\4 
CL 1 (1 J 

26 IF(1— Ai\Gf l"‘"7)27*'^C*tr' ^ 

27 t.'NNU = ANNU«G INC 
uC 1 CJ J 

LN’r4U = AN\U »S I N7 
GC TO 3 
29 l-NNU=Ah4,\U 

J. cNRt=W Au * ( 1 . r2)/(( . A02S6FTF ) t32.b3? > 

IF(LOF ) .21,2? 

21 TFF=TT-( .8 »Xl.wTC) 

GC 1" 23 

22 TFF-TT-( .Ov «XLN’TO) 

23 CNRr:= ( ft’ Alj < ( 1 . F ^: ))/((' . A 6 2 '8 * f F f ) + 3 ' . 8 8 . ' ) 

FHO= ( ( P 1 t 1 3 . r >e*Pfl ) * ' . 1 f- I 5 ) / ( f 3 . ( 1 FF I '»59 . C 9 ) ) 

FR I F AC- ( 1 . 3 8A »RHL.«Pf * ( 1 . 'f 8 ))/({' AiJ/ . 1 6 7 3 > * >^ V ) 

I F(l._ GF > . .31.32 

J 1 I r ( L A F ) 2 , 3 3 , 3 A 

32 IF(LAF)^ .3‘:-.36 

33 ARGrAC=(4. 74 « ( 1^4Rt **2 ) ) / ( I . : ^ 1 ) - ( 1 . o 74 n ( h,N' • > ) / ( 1 . 8 ) F 1 . 1 6 ':4 ft 

331 I'lAFAC^ .88 3 

cC IF 37 

34 A 4GF LC= ( 3 . 74 7* < FMkf c r 2 ) ) / ( 1 . ''1 > - ( 1 . 7 74 » ( .3 \XE ) ) / ( 1 . '8)-H.l'F.l8 

*341 P1AFAC= .8278 

GC Ti; 3/ 

35 ARGFAC= .94 f ( 3 . 4 j ' « 8NRF ) / 1 , ’ o- ( . I 2 9 * ( 9NR'' » * 2 ) ) / 1 . , t 9 

351 CIArAC= .883. 

C l; r ( i 3 7 

36 ARGF AC- ( 4 . 6 3* ( *«))/!. I 1- ( ■>. 36 ;* ) / 1 . 'F 4 I . r.6 

361 CIAFAC= .327” 

37 TNNl = ( (3 j 3!' = « * .6 )*,. 92'LARGiAC»(((C.r'.»3. «( .16*Tr+ll?.8)>/((l.'LA 

1»»( . ' 4. . 1344)))** . 33.3) 

FTHi FA- •NNU/TNr4b 

36 CNMU=(TNfsU* . 7> /7* ( ■ NRt. » « .3 ) * 0 I Al- A C ) / A < GF A C 

XR I I h ( 6 . 4 ) IRCN . fnnU , TNNU . C MNU. r FH'. FA . :.,NRr ,FR I FAC , CNR., 
f. T = - T- I 
2 1 CCl'iriNU 

crop 

ENH 



I 



/ 
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APPENDIX E 
Experimental Data 

This Appendix presents a listing of all test data. The 
data is listed sequentially by run code number^. No lobels 
are proviaed on the data listings, accordingly the follov/ing 
information is provided: 



Column 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 



Description 

Run code number for data point identification, 
four digits, (RNCD), 

Pressure drop across the thin plate orifices, 
inches of oil, three or four digits (DHF) 
Pressure ahead of the thin plate orifice, inches 
of water, three or four digits (Pj^) 

Barometric pressure, inches of mercury, four 
digits (PB) 

Pressure drop across the tube bank inches of 
Mercury, four digits (DP) 

Inlet temperature, degrees Farenheit, four 
digits (TI) 

Exhaust temperature, degrees Farenheit, four 
digits (TE) 

Tube surface temperature, degrees Farenheit, 
four digits, (TT), 

Orifice diameter, incties (IDO) 

Specific humidity, grains per pound (SPHU) 
Geometry indicator 'LGF) 

Angle indicator (LANGLE) 

Arrangement indicator (LAF) 



1401 


1,40 


1.00 


30.46 


0.100 


1402 


1.50 


1.40 


29 . 70 


0.340 


1403 


1.70 


1.20 


30.46 


0.120 


1404 


1.70 


i.ro 


29.70 


0.100 


1405 


2.00 


1.20 


29.95 


0.160 


1406 


2.00 


1.60 


30.46 


0.130 


1407 


2.50 


1.20 


30.46 


0.120 


1408 


2.50 


1.80 


29.70 


0.490 


1409 


2.70 


1.60 


29.95 


0.130 


1410 


2.80 


1.90 


29.70 


0.160 


1411 


2.90 


2.00 


30.46 


0.130 


1412 


3.30 


2.00 


29.95 


0.160 


1413 


3.90 


2.40 


29 . 70 


0.210 


1414 


4.10 


2.50 


29.95 


0.210 


1415 


4.20 


3.00 


30.46 


0.210 


1416 


4 . 3 u 


3.20 


29.70 


0.800 


1417 


4.90 


3.00 


29 . 95 


0.230 


1418 


4.90 


3.00 


30.46 


0.230 


1419 


5.70 


3.60 


29.95 


0.310 


1420 


6.20 


4.00 


29 . 70 


0.360 


1421 


6.50 


4.00 


30.46 


0.570 


1422 


7.80 


5.00 


29.70 


0.390 


1423 


8.30 


6.20 


29.70 


1.550 


1424 


8.30 


5.20 


30.46 


0.390 


1425 


8.50 


6.40 


29.70 


1.630 


1426 


9.60 


6.50 


30.46 


0.570 


1427 


11.40 


7.60 


29.70 


0.620 


1426 


13.10 


9.20 


30.46 


0.670 


1429 


16.00 


9.60 


29.70 


0.830 


1430 


16.10 


10.50 


30.46 


0.830 


1431 


18.40 


12.20 


30.46 


0.830 


1432 


20.50 


13.60 


30.46 


1.040 


14 33 


24.00 


16.00 


30.46 


1.220 


1434 


25 . 5 G 


16.80 


29.70 


1.300 


1435 


29.10 


19.20 


30.46 


1.400 


1436 


3 3.00 


22 . CO 


30.46 


1.610 


1601 


1.10 


0.50 


30.48 


0.200 


1602 


1.60 


1.20 


30.48 


0.400 


1603 


2.60 


1.60 


30.48 


0.500 


1604 


2.90 


2.40 


29.82 


1.100 


1605 


3 . 3 C 


2.00 


30.48 


o.soo 


1606 


4.10 


2.60 


30.48 


0.900 


1607 


5.00 


3.30 


30.48 


1.200 


1608 


6.00 


4.00 


30.48 


1.400 


1609 


6.40 


4.20 


29.82 


1.900 


1610 


7 . CO 


4 . 60 


3 Q .48 


1.700 


1611 


7.20 


5.00 


29.82 


2.100 


1612 


8.10 


5.40 


29.82 


2.400 


1613 


8.30 


5.40 


30.48 


1.700 


1614 


9 . 4 C 


6.40 


3 C .48 


2.500 


1615 


10.00 


6.60 


29.82 


2.600 


1616 


11.60 


7.40 


29.82 


3.2 00 


1617 


11.60 


6.80 


30.48 


2.900 


1618 


12.60 


8.60 


30.48 


4.000 


1619 


12 . 9 C 


8.80 


29.82 


3.600 



.0 


176.0 


212.0 


3 


77.0 


0 


4 


0 


.0 


166.0 


212.0 


4 


87,0 


0 


4 


0 


.0 


174 . C 


212.0 


3 


77.0 


0 


4 


0 


.0 


176.0 


212.0 


3 


94.0 


0 


4 


0 


.0 


165.0 


212.0 


3 


70.0 


0 


4 


0 


.0 


163.0 


212.0 


3 


77.0 


0 


4 


0 


.0 


174.0 


212.0 


3 


77.0 


0 


4 


0 


.0 


162.0 


212.0 


4 


87.0 


0 


4 


0 


.0 


163.0 


212.0 


3 


70.0 


0 


4 


0 


.0 


167.0 


212.0 


3 


94.0 


0 


4 


0 


.0 


163.0 


212.0 


3 


77.0 


0 


4 


0 


.0 


164.0 


212.0 


3 


70.0 


0 


4 


0 


.0 


176.0 


212.0 


3 


94.0 


0 


4 


0 


.0 


163.0 


212.0 


3 


70.0 


0 


4 


0 


.0 


163.0 


212.0 


3 


77,0 


n 


4 


0 


.0 


162.0 


212.0 


4 


87.0 


0 


4 


0 


.0 


162,0 


212.0 


3 


70.0 


0 


4 


n 


.0 


162.0 


212.0 


3 


77.0 


0 


4 


0 


.0 


165.0 


212.0 


3 


70.0 


0 


4 


0 


.0 


168.0 


212.0 


3 


94 .0 


0 


4 


0 


.0 


' 160.0 


212.0 


3 


77.0 


0 


4 


0 


.0 


167,0 


212 .^/ 


3 


94.0 


0 


4 


0 


.0 


155.0 


212.0 


4 


87.0 


0 


4 


0 


.0 


157,0 


212.0 


3 


77.0 


0 


4 


0 


.0 


157,0 


212.0 


4 


87.0 


0 


4 


0 


.0 


160,0 


212.0 


3 


77.0 


0 


4 


0 


.0 


162.0 


212.0 




94.0 


0 


4 


0 


.0 


161.0 


212.0 


3 


77.0 


0 


4 


0 


.0 


150.0 


212.0 


3 


94.0 


0 


4 


0 


.0 


159.0 


212.0 


3 


77.0 


0 


4 


0 


. 0 


158.0 


212.0 


3 


77.0 


0 


4 


0 


• 0 


159.0 


212.0 


3 


77.0 


0 


4 


0 


.0 


136.0 


212.0 


3 


77.0 


0 


4 


0 


.0 


156.0 


212.0 


3 


94,0 


0 


4 


0 


.0 


156.0 


212.0 


3 


77.0 


0 


4 


0 


.0 


155.0 


212.0 


3 


77.0 


0 


4 


0 


.0 


171.0 


212.0 


3 


62,0 


0 


6 


0 


.0 


164.0 


212.0 


3 


62,0 


0 


6 


0 


.0 


165.0 


212.0 


3 


62.0 


0 


6 


0 


.0 


166,0 


212.0 


3 


46.0 


0 


6 


0 


.0 


165.0 


212.0 


3 


62.0 


0 


6 


0 


.0 


164.0 


212.0 


3 


62.0 


0 


6 


0 


.0 


161.0 


212.0 


3 


62,0 


0 


6 


0 


.0 


160.0 


212.0 


3 


62.0 


0 


6 


0 


.0 


162.0 


212.0 


3 


46.0 


0 


6 


0 


.0 


156.0 


212.0 


3 


62,0 


0 


6 


0 


.0 


161.0 


212.0 


3 


46.0 


0 


6 


0 


.0 


161 .0 


212 . C 


3 


46.0 


0 


6 


0 


.0 


158.0 


212.0 


3 


62.0 


0 


6 


0 


.0 


155.0 


212.0 


3 


62,0 


0 


6 


0 


.0 


160.0 


212.0 


3 


46.0 


0 


6 


0 


.0 


160.0 


212.0 


3 


46 . 0 


0 


6 


0 


.0 


156.0 


212.0 


3 


62.0 


0 


6 


0 


.0 


154.0 


212.0 


3 


62.0 


0 


6 


n 


.0 


159,0 


212.0 


3 


46.0 


0 


6 


n 



80 

70 

78 

79 

72 

72 

78 

72 

72 

79 

72 

72 

79 

73 

73 

72 

73 

74 

74 

79 

74 

79 

78 

75 

80 

78 

80 

So 

82 

82 

83 

64 

65 

84 

87 

63 

74 

75 

76 

61 

76 

77 

78 

78 

62 

78 

81 

82 

78 

80 

82 

83 

81 

83 

63 



I 




1620 


14.60 


9.80 


29.82 


4.200 


1621 


15.20 


10.60 


30.48 


4.000 


1622 


16.50 


11.00 


29.82 


4.5 00 


1623 


17.90 


12.30 


29 . 82 


5.000 


1624 


18 . 6 C 


12.60 


29.82 


5.000 


1625 


18.60 


12.80 


30.48 


5.300 


1626 


20.50 


14.00 


29.82 


5.700 


1627 


22.30 


15.20 


29.82 


5.900 


1626 


23.40 


15.90 


30.48 


6 . ICO 


1629 


24.40 


16.60 


29.82 


6.500 


1630 


26.00 


17.80 


30.48 


7.000 


1631 


26.40 


18.20 


29.62 


7.000 


1701 


1.20 


0.80 


29.61 


0.400 


1702 


2.20 


1.60 


29.61 


0 . 60 C 


1703 


3.30 


2 . O ' 0 


29.61 


1.200 


1704 


4.20 


2.60 


29.61 


1.500 


1705 


5.30 


3.50 


29.61 


2.000 


1706 


6.30 


4.10 


29.61 


2.4 00 


1707 


7.20 


4.80 


29.61 


2.500 


1708 


8.40 


5.60 


29.61 


2.900 


1709 


9.30 


6.30 


29.61 


3.400 


1710 


10.70 


8.00 


29.61 


3.900 


1711 


12.50 


8.90 


29.61 


4.300 


1712 


14.40 


10.10 


29.61 


4.800 


1713 


16.40 


11.60 


29.61 


5.500 


1714 


18.20 


13.00 


29.61 


6.000 


1715 


20.10 


14.30 


29.61 


6 . v 800 


1716 


23.20 


16.50 


29.61 


7,800 


1717 


26.70 


18.50 


29.61 


8.900 


1718 


2 8.90 


19.90 


29.61 


9.700 


1901 


2.80 


2 . CO 


30.34 


0.310 


19 C 2 


3.90 


3.90 


30.34 


1.740 


1903 


6.60 


6.70 


30.34 


2 . S 50 


1904 


9.00 


6.13 


30.34 


3.000 


1905 


10.90 


10.6 0 


30.34 


4.540 


1906 


16 . GO 


11.70 


30.34 


4.000 


1907 


27.20 


19.40 


30.34 


5.220 


1901 


25.00 


29.40 


30.08 


1.000 


1902 


20.5 


24.40 


30.08 


1.00 


1903 


18.0 


22.6 


30.08 


1.0 


1904 


9.0 


12.2 


30.08 


1.0 


1905 


20.05 


22.0 


29.75 


1.0 


1906 


17.25 


19.12 


29.75 


1.0 


1907 


14 . 3 B 


17.4 


29.75 


1.0 


1908 


7.88 


10.42 


29.75 


1.0 


1909 


6.25 


5 . 2 ? 


29.75 


1.0 


1910 


5.5 


8.7 


29 . 75 


1.0 


1911 


3.66 


0.0 


29.75 


1.0 


1912 


20.62 


25.0 


29.68 


8.78 


1913 


17.12 


21.0 


29.68 


7,37 


1914 


12.63 


15.5 


29.68 


6.25 


1915 


11 . 1 ? 


13.5 


29.68 


5.12 


1916 


5.5 


6.87 


29.68 


2.63 


1917 


4 . 7 :> 


6.2 5 


29.68 


2.37 


1918 


1.63 


2.50 


29.68 


0.88 



.0 


158.0 


212.0 


3 


46.0 


0 


6 


0 


.0 


154.0 


212.0 


3 


46.0 


0 


6 


0 


.0 


156.0 


212.0 


3 


46.0 


0 


6 


0 


.0 


156.0 


212.0 


3 


46.0 


0 


6 


0 


.0 


157.0 


212 .^ 


3 


46.0 


0 


6 


0 


.0 


153.0 


212.0 


Q 


46.0 


0 


6 


0 


,0 


155.0 


212.0 


3 


46.0 


0 


6 


0 


,0 


157.0 


212.0 


3 


46.0 


0 


6 


0 


.0 


154.0 


212.0 


3 


62.0 


0 


6 


0 


.0 


155.0 


212.0 


3 


46.0 


0 


6 


0 


.0 


154.0 


212.0 


3 


62.0 


0 


6 


0 


,0 


154.0 


212.0 


3 


46.0 


0 


6 


0 


.0 


166.0 


212.0 


3 


42.0 


0 


7 


0 


.0 


160,0 


212.0 


3 


42.0 


3 


7 


n 


.0 


159.0 


212.0 


3 


42.0 


0 


7 


0 


.0 


159.0 


212.0 


3 


42.0 


0 


7 


0 


.0 


158,0 


212.0 


3 


42.0 


0 


7 


0 


.0 


156.0 


212.0 


3 


42.0 


0 


7 


0 


. r ' 


157,0 


212.0 


3 


42.0 


0 


7 


0 


.0 


158,0 


212.0 


3 


42.0 


0 


7 


0 


.0 


156.0 


212.0 


3 


42.0 


0 


7 


0 


,0 


156,0 


212.0 


3 


42.0 


0 


7 


0 


. 0 


154,0 


212.0 


3 


42.0 


0 


7 


0 


,0 


153,0 


212.0 


3 


42.0 


0 


7 


0 


. 0 


153.0 


212.0 


3 


42.0 


0 


7 


n 


.0 


153.0 


212.0 


3 


42.0 


0 


7 


0 


.0 


153.0 


212.0 


3 


42.0 


0 


7 


0 


,0 


15 3.0 


212.0 


3 


42.0 


0 


7 


0 


.0 


153.0 


212.0 


3 


42.0 


0 


7 


0 


,0 


15 3.0 


212.0 


3 


42.0 


0 


7 


n 


.0 


169.0 


212.0 


3 


84.0 


0 


9 


0 


.0 


156.0 


212.0 


4 


84,0 


0 


9 


0 


.0 


152.0 


212.0 


4 


84.0 


0 


9 


0 


.0 


157.0 


212.0 


3 


84.0 


0 


9 


0 


,0 


15 2,0 


215.0 


4 


04.0 


0 


9 


0 


.0 


15 5.0 


210.0 


3 


84,0 


c 


9 


0 


,0 


154,0 


212.0 


3 


84.0 


0 


9 


0 


.0 


151,0 


204.0 


4 


50.0 


0 


9 


0 


.0 


145.0 


208.0 


4 


50 . r 


0 


Q 


0 


.0 


145 .0 


210 . C 


4 


50 . 


0 


9 


o 


.0 


144,0 


210.0 


4 


50.0 


0 


9 


0 


.0 


147.0 


215.0 


4 


68.0 


0 


9 


0 


.0 


146,0 


210.0 


4 


68.0 


0 


9 


0 


.0 


146,0 


205.0 


4 


68.0 


0 


9 


0 


,0 


149,0 


210.0 


4 


68.0 


0 


9 


0 


.0 


150,0 


212.0 


4 


68.0 


0 


9 


0 


.0 


151.0 
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APPENDIX F 

Reduced Experimental Data 

\ 

This Appendix consists of a listing of the results of 
the data reduction. The column headings are self explanatory. 
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APPENDIX G 

Convection Heat Transfer Regression Analysis 
A regression analysis is applied to the results of the 
data reducing computer program, the purpose of v;hich is to 
fit a curve of the generalized type, Y = A + BX^ through 
the reduced data. The independent variable X is the exper- 
imental value of the Reynolas Number through the test section, 
defined throughout this text as BNRE. Tne indepenaent 
variable Y is the ratio of the Nusselt Number as measured 
experimentally to the Nusselt Number computed from the 
modified Grimison equation, defined by FTHETA, The results 
of this regression analysis will be the constants A and B, 
and the exponent K resulting in an equation relating FTHETA 
to Reynolds Number BNRE, within the range of the experimental 
data. From this equation, particularized for tne angle of 
inclination and the tube geometry, the angle factor FTHETA 
may be computed, or taken from design-oriented curves to apply 
as an additional factor to the modified Grimison Equation, 



The analysis was carried out using a prepared program 
available at the Applied Mathematics Laboratory of the David 
Taylor Model Basin, The program is written in Bell FORTRAN II 
for the IBM 7090 computer. Running time is 2,5 minutes. The 
program is documented in DTiMB Report 20 3 7, 

A brief description of the program is included in this 
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A maximum of 35 data points may be entered as input for 
any one case. The program will fit a function of the form 



squares. The line of regression will be calculated and plotted 
three times using three sets of data points, the first being 
the entire original set of data points, the second being 
the set of points within two standard errors, ano tne third 
using those points which fall within one standard error. The 
standard error is defined as: 



where: Y average is the average ordinate of all the points. 

Yi is the general ordinate when i is indexed from 1 to N. 

N is the number of points. 

The method used by the program is the method of least 
squares. An approximation to the exponent M is made by using 
the first three and the last three points, averaging the 
ordinates and abcissas of these points to derive an average 
point, X average. Four points, (two greater than and two less 
than) closest to X average are taken. From these four points, 
a new average point is derived. The slopes of the upper and 
lower portion of the first approximate line of regression are 
then determined. Dividing the lower slope by the upper slope, 
an expression exists from which the exponent M may be approximated. 



Y = A+B X 



M 



to the data points entered by the method of least 
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Xs = X^.X2.X3 



Ys = Y^.Y2.Y3 



X- = X +X o 

L n n-1 n-2 



Y^ = Y +Y t+Y „ 
L n n-1 n-2 



X avg , 






The four values of X closest to X avg. are located, 



X„ = X. „ +X. T + X. + X. + 1 
C 1-2 1-1 i 1 



Y„ = Y. ^ + Y. T +Y. +Y. T 
C 1-2 1-1 1 1+1 



in which X^ = X avg. 



Then the exponent M may be determined in 
the following manner: 

M 

Y = A + BX 

M-1 

= BFiX 



dx 
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Two slopes may be written from Fig. A 



(1) -Y3 



= BM 






2 / 

where is the slope at + Y^ 



\ 



""c - ""s 



(2) Y^ - Y^ 
Xl - X^ 



= BM 






M-1 



where 





the 



slope 



at X^ + X^ 
2 



Dividing equation (2) by equation (1). 



(Xl - X^) (Y^' - Y3) 





1 



Now let (Y^ “ 



Let . X^ ^ ^ 

Xc -H X3 

D = E 

V 

LnD = (M-1) In E 

Then M — In D 
In E 

which provides an approximation to M. 

The initial approximation of the exponent is testing for 
accuracy by considering seven valut-S of the expon nt, M being 
the center value and the others 1/4 M from this central value in 
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increasing ana decreasing increments. 

For each value of the exponent, the values of A, ana B 
are computed, and the sums of the squares of the derivations 
are minimized. The exponent, M, is then chosen from the set 
of A, B, and M that has yielded the rninir.ium sum of the deriva- 
tions squared, 

M 

The normal equations for the function Y = A + BX v^ith 
A, B, and M as the regression parameters. 



These equations are solved for A ana B by Cramer’s Rule. 
The final test for the proper exponent is maae by examining 
the coefficients of the three normal equations, i.e., three 
equations in two unknowns. Since A and B are unique, two of 
the three equations must be linearly dependent ana the 
determinant of the coefficients will be identically equal to 
zero. 
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This is true because the coefficients of one equ tion 
must be equal or proportional to the coefficients of another. 

It can now be seen that successive iterations basea on estimates 
of M until the determinants of the coefficients are equal to 
zero to finalize the value of the coefficient M, 

The exact values of A and B can then be calculated using 
Cramer’s rule. 

The output of the program appears in two forms. One form 
is the regular output of the 7090 program which states the 
A and B coefficients and the exponent M for each of the three 
lines of regression. The output of a typical analysis is 
shown on pagei 3 J|for examination. Note that with the lines of 
regression for points within two and one stanuara error, those 
points which were net used in the determination of the line 
are listed with the results. This list of points is cumulative 
through tne two cases. The identifying information v;hich 
accompanies each uata point is placed in the computer along 
with the data at the choice of the user of the program. The 
second form of output is through the SC 4020 Charactron, The 
coaing for this output is contained in a subroutine package 
named AMPFOS and the form of the output is microfilm copies 
of the actual lines of regression. Although these plots 
have little value from the practical standpoint, due to the 
selection by the program's coding of rather cumbersome scales, 
the plots do indicate the validity of the curve fit ano permit 
selection of one of the three lines as being mo^t representative 
of the points. This output form is shown on page 1^5. 
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Results of Data Reduction 

Figures 27 thru 60 which follow ure plots of Nusselt 
Number versus Reynolas Number for individual geometries and 
angles. Curves for 0 = 90, 75, 60, 45 and thirty degrees are 
plotted using experimentally derived data, other curves for 
each geometry are also plotted using the Colburn ane Grimison 
equations. 
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PUNS 1^01- 1^36 




FIG 2 9 
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nUNS 1601-1631 




FIG 20 
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FIG. 43 
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FIG. 44 
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FIG. 50 



161 



APPENDIX H 

' Analysis of Frictional Resistance Results 

Figures 51 thru 66 are plots of the computed friction 

factor versus Reynolds Number for individual geometries and 

/ 

angles of inclination. The uncertainty of data for Reynolds 
Numbers lower than 2500 was such that no points were used in 
the analysis. The process of curve fitting used was the 
method of least squar^. 

Least Squares Analysis 

The frictional resistance data is fittea to an equation 
of the form: 

F YPO 

FRIFAC = CONST (CNRE) 

using the method of least squares. The logarithmic equivalent 
of the equation is taken: 

log FRIFAC = log CONST + EXPO log CNRE 
and the normal equations developed by the usual method of 
multiplying by the coefficients of the variables CONST and EXPO 
in each case for each data point and summing. The normal 
equations so developed are; 

log FRIFAC = (log CNRE) EXPO + THUMB log CONST 

(log FRIFAC log CNRE) = (log CNRE)^ EXPO + (log CNRE)log CONST 
These normal equations can be solved simultaneously for the log 
of CONST and EXPO: 

EXPO = (THUMB) (log FRIFAC) (log CNRE)-(loq CNRE) (log FRIFAC) 

(THUMB) (log CNRE)^ + (log CNRE)^ 

XCON = (log FRIFAC) - (log CNRE) (EXPO) 

(THUMB) ~ 

Then: CI3NST = antilog XCON 

The parameters necessary for expression of friction factor as a 
function of Reynolas Numbers are thus developed. 
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RUNS 2601-2626 
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Frictional Resistance Analysis Computer Program Description 

The computer program described herein is a FORTRAN II 
program compatible with trie IBM 1620 data processing system. 

The function of the program is to fit a curve whereby friction 
factor is expressed as a power of Reynolds Number using the 
method of least squares. While variable names, input, output 
and format statements are particularized to the problem at hana 
the methoa is completely general and can be used in any 
application where a logarithmic least squares analysis is 
desired with only minor program moaif ications . It would be 
possible to include tnis program in the basic Uata reduction 
program to permit reduction ana analysis in one computation. 

Overall control of the program i s by a DO loop inaexed by 
the number of sets of oata being analyzed. For each uata set 
another DO loop inaexed by the number of data points in the 
data set controls the calculations. For each data set the 
progrcim first sets the summations equal to NUME , then accepts 
as input the friction factor and Reynolds number for the first 
data point. The natural logarithms of the input quantities, 
the product of the two logarithms ana the square of the 
logarithm of Reynolds Number are coraputea. The summations 
necessary to determine the coefficients of tne normal equations 
are computed with variable names being changed as necessary 
for proper program operation. The program tiien loops to repeat 
the calculations described above for each data point until all 
data points in the set have be<- n calculated. Tne values of 
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the exponent, hXPO , ana the natural logarithm of the coefficient, 
XCON can then be computed for the particular data set. As no 
antilogarithm function is available in the FORTRAN II library 
the value of the coefficient, CONST, is calculated by: 

CONST = 

When tiie values of the coefficient and exponent have been 
determined an output consisting of these quantities plus an 
identifying coae , and the logarithm of the coefficient is 
printed . 



The quantities which have been computed so far represent the 
parameters necessary for expressing friction factor as a function 
of Reynolus Number for a particular geometry ana angle of 
inclination. In order to correlate the values calculated with 
the equation aeveloped using these parameters with test results, 
values of Reynolus Number from three thousand to eight thousand 
are inserted into the equation and the corresponaing friction 
factor computed. The results are printed. 

VJhen the calculation of the friction factors is completed 
the program loops to perform the calculations for trie next data 
set until the limiting index is reached. 



A program block diagram, listing, sample output, and output 
summary follow. 

Computations using this program were made at C.W. Post 
College; data was prepared at Webb Institute, 
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C FORTRAN TWO PROGRAM TO CORRELATE FRICTIONAL RESISTANCE DATA BY 

C LEAST SUCARES 

1 FORKAT(U) 

2 FORMAT! 2 14) 

3 FORMAT ( F9.6 , 2X , P9. 2 ) 

4 FORMAT ( Ihl , 2X, 57HRUN CODE EXPONENT CONSTANT LOGlCUNo 

P ITANT) ) 

b FORMAT ( 1H0,4X, I2,7X,F10.6,5X,F10.6,5X,F10.6) 

9 FORMAT (1H0,6X,F8.2, 13X,F9.6) 

10 FORMAT { 1H0,2X, 37hREYN0LDS NUMBER FRICTION FACTOR ) 

READ l,KASeS 

DO 7 MA=1,KASES 

SUMFLN=0.0 

SUMREL=O.U 

SMSRG=0.0 

SMFGRE=0.0 

READ 2,NUMB,RNCD 

DO 6 NU=1,NUMB 

READ 3,FRIFAC,CNRE 

FLG=LOGF(FRIFAC) 

RFLG=LOGF (CNRE ) 

FLGREL=FLG*RELG 

SRELG=RELGY*2 

sumflg=sumfln+flg 

SUMFLN=SUMFLG 

smrelg=slmrel+relg 

sumrel=smkflg 

SMSREL=SMSRG+SRELG 
SMSRG=SMSREL 
SMFLRL=SMFGRE+FLGREL 
SMf-GRE = SMFLRL 

6 CONTINUE 
THUMB=NUMB 

EXP0=( THLMB*SMFGRE-SUMREL=-SUMFLN ) / ( ThCMB'-^SMSRG- ( SLMR E L 2 ) ) 

XCCN=( ( SUMFLN-! SUMREL^EXPO ) ) /THUMB) 

CPNST = 2.7ia28-=<'XC0N 
PRINT 4 

PRINT 5, RNCO, EXPO, CONST, XCUN 
PRINT 10 
DO 8 NRE=3,8 
DNRE=NRE 

ENRE=DNKE*( 1 ,0E3) 

F R I C T =CONS T--:‘ ( ENR E -y =:^E XPO ) 

PRINT 9,ENRE,FRICT 
8 CONTINUE 

7 CONTINUE 
STOP 
END 
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kUN CODE EXPONENT 


CONSTANT 


L00( CONSTANT ) 


43 .060958 


.033258 


-3.403A55 ■ 


REYNOLDS NUMBER 


FRICTION factor 




3000.00 


.05^187 




4000.00 


.0551^6 




6000.00 


. 055901 




6000.00 


.056526 




7000.00 


.057059 




8000.00 


.057526 





STOP 






16 

17 

19 

24 

26 

29 

33 

34 

36 

39 

43 

44 

45 

49 
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Summary of Frictional Resistance Analysis 



Exponent 

- 0.043053 

0.196590 

- 0.045756 

- 0.079225 

- 0.024007 

- 0.119514 

- 0.081586 

0.057896 

- 0.089713 

- 0.028274 

- 0.086865 

0.060968 

0.210122 

0.053331 

- 0.055798 



Constant 

0.065543 

0.045465 

0.445318 

0.214984 

0.100758 

0.037337 

0.261858 

0.040784 

0,270171 

0.1877P7 

O.U31623 

0.033258 

0.017001 

0.074656 

0.297399 



j 
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Run Code 


14 


16 


17 


19 


Reynolds 

Number 


Friction 

Factor 


Friction 

Factor 


Friction 

Factor 


Friction 

Factor 


3000 


0.046433 


0.219407 


0.308722 


0.114006 


4000 


0.045861 


0.232174 


0.304685 


0 J11437 


5000 


0.045422 


0.242585 


0.301590 


0.109484 


6000 


0.045067 


0.251438 


0.299084 


0.107914 


7000 


C. 044769 


0.259174 


0.296982 


0.106604 


8000 


0.044513 


0.266068 


0.295173 


0.105482 


Run Code 


24 


26 


29 




Reynolus 

Number 


Friction 

Factor 


Friction 

Factor 


Friction 

Factor 




3000 


0.083138 


0.097209 


0.136263 




4000 


0.082566 


0.110610 


0.133102 




5000 


0.082125 


0.103329 


0.130701 




6000 


0.081766 


0.105605 


0.128771 




7000 


0.081464 


0.107569 


0.127162 




8000 


0.081204 


0.109299 


0.125784 




Run Code 


33 


34 


36 


39 


Reynolds 

Number 


Friction 

Factor 


Friction 

Factor 


Friction 

Factor 


Friction 

Factor 


3000 


0.064834 


0.131732 


0.149752 


0.215310 


4000 


0.065923 


0.128376 


0.148539 


0.209996 


5000 


0.066780 


0.125831 


0.147605 


0.205965 


6000 


0.067489 


0.123790 


0.146846 


0.202728 


7000 


0.068094 


0.122090 


0.146207 


0.200032 


8000 


0.068622 


0.120636 


0.145656 


0.197725 
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Run Code 


43 


44 


46 


49 


Reynolds 

Number 


Fricti on 
Factor 


Fricti on 
Factor 


Friction 

Factor 


Friction 

Factor 


3000 


0.054187 


0.091434 


0.114421 


0.190248 


4000 


0.055145 


0.097132 


0.116190 


0.187219 


5000 


0.055901 


0.101795 


0.117581 


0.184902 


6000 


0.056526 


0.105770 


0.118730 


0.183031 


7000 


0.051059 


0.109252 


0.119710 


0.181463 


8000 


0.057526 


0.112361 


0.120568 


0.180116 



Friction Factor (f) 
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